Dynamical Casimir effect for gravitons in bouncing braneworlds 
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We consider a two-brane system in five-dimensional anti-de Sitter space-time. We study particle 
creation due to the motion of the physical brane which first approaches the second static brane 
(contraction) and then recedes from it (expansion). The spectrum and the energy density of the 
generated gravitons are calculated. We show that the massless gravitons have a blue spectrum 
and that their energy density satisfies the nucleosynthesis bound with very mild constraints on the 
parameters. We also show that the Kaluza-Klein modes cannot provide the dark matter in an anti- 
de-Sitter braneworld. However, for natural choices of parameters, backreaction from the Kaluza- 
Klein gravitons may well become important. The main findings of this work have been published 
in form of a Letter [R.Durrer and M.Ruser, Phys. Rev. Lett. 99, 071601 (2007), arXiv:0704.0756]. 
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I. INTRODUCTION 

In recent times, the possibility that our observed 
Universe might represent a hypersurface in a higher- 
dimensional space-time has received considerable 
attention. The main motivation for this idea is the 
fact, that string theory |l], 0], which is consistent only 
in ten spac-etime dimensions (or 11 for M-theory) 
allows for solutions where the standard model particles 
(like fermions and gauge bosons) are confined to some 
hypersurface, called the brane, and only the graviton can 
propagate in the whole space-time, the bulk |l| . Since 
gravity is not well constrained at small distances, the 
dimensions normal to the brane, the extra dimensions, 
can be as large as 0.1mm. 

Based on this feature, Arkani-Hamed, Dimopoulos 
and Dvali (ADD) proposed a braneworld model where 
the presence of two or more flat extra-dimensions can 
provide a solution to the hierarchy problem, the problem 
of the huge difference between the Planck scale and the 
electroweak scale ||, || . 

In 1999 Randall and Sundrum (RS) introduced a model 
with one extra dimension, where the bulk is a slice of 
five-dimensional anti de-Sitter (AdS) space. Such curved 
extra dimensions are also referred to as warped extra 
dimensions. While in the RS I model Q with two flat 
branes of opposite tension at the edges of the bulk the 
warping leads to an interesting solution of the hierarchy 
problem, it localizes four-dimensional gravity on a single 
positive tension brane in the RS II model M. 
Within the context of warped braneworlds, cosmological 
evolution, i.e., the expansion of the Universe, can be 
understood as the motion of the brane representing our 
Universe through the AdS bulk. Thereby the Lanczos- 
Sen-Darmois-Israel-junction conditions [g, ||, [l(| [Til , 
relate the energy-momentum tensor on the brane to the 
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extrinsic curvature and hence to the brane motion which 
is described by a modified Friedmann equation. At low 
energy, however, the usual Friedmann equations for the 
expansion of the Universe are recovered [ 12|, [l3| . 

Since gravity probes the extra dimension, gravita- 
tional perturbations on the brane, i.e. in our Universe, 
carry five-dimensional effects in form of massive four- 
dimensional gravitons, the so-called Kaluza-Klein (KK) 
tower. Depending on the particular brane trajectory, 
these perturbations may be significantly amplified lead- 
ing to observable consequences, for example, a stochastic 
gravitational wave background. (For a review of stochas- 
tic gravitational waves see 63] •) This amplification 
mechanism is identical to the dynamical Casimir effect 
for the electromagnetic field in cavities with dynamical 
walls (moving mirrors); see [|l5| |l7j and references 
therein. In the quantum field theoretical language, 
such an amplification corresponds to the creation of 
particles out of vacuum fluctuations. Hence, in the same 
way a moving mirror leads to production of photons, 
the brane moving through the bulk causes creation of 
gravitons. Thereby, not only the usual four-dimensional 
graviton might be produced, but also gravitons of the 
KK tower can be excited. Those massive gravitons are of 
particular interest, since their energy density could dom- 
inate the energy density of the Universe and spoil the 
phenomenology if their production is sufficiently copious. 

The evolution of cosmological perturbations under 
the influence of a moving brane has been the subject of 
many studies during recent years. Since one has to deal 
with partial differential equations and time-dependent 
boundary conditions, the investigation of the evolution 
of perturbations in the background of a moving brane is 
quite complicated. Analytical progress has been made 
based on approximations like the "near brane limit" and 
a slowly moving brane @ |(| §l| . 
The case of de Sitter or quasi-de Sitter inflation 
on the brane has been investigated analytically in 
|E1 51 E2I pi E6|. In pi it is demonstrated that dur- 



2 




FIG. 1: Two branes in an AdSs spacetime, with y denot- 
ing the fifth dimension and L the AdS curvature scale. The 
physical brane is on the left at time dependent position yb(t). 
While it is approaching the static brane its scale factor is 
decreasing and when i t m oves away from the static brane it 
is expanding [cf. Eq. (2.3)]. The value of the scale factor of 



the brane metric as function of the extra dimension y is also 
indicated. 



ing slow-roll inflation (modeled as a period of quasi-de 
Sitter expansion) the standard four-dimensional result 
for the amplitude of perturbations is recovered at low 
energies while it is enhanced at high energies. 
However, most of_the_ effort has gone into numerical 
simulations 



35, 36fl, in 



2% |2j, |29|, £0|, y £3, |33|, 

particular in order to investigate the high-energy regime. 
Thereby different coordinate systems have been used 
for which the brane is at rest, and different numerical 
evolution schemes have been employed in order to solve 
the partial differential equation. 

In this work we chose a different way of looking at 
the problem. We shall apply a formalism used to 
describe the dynamical Casimir effect to study the 
production of gravitons in braneworld cosmology. This 
approach and its numerical implementation offers many 
advantages. The most important one is the fact that 
this approach deals directly with the appearing mode 
couplings by means of coupling matrices. (In |]l9| a 
similar approach involving coupling matrices has been 
used. However, perturbatively only, and not in the 
complexity presented here.) Hence, the interaction 
between the four-dimensional graviton and the KK 
modes is not hidden within a numerical simulation but 
can directly be investigated making it possible to reveal 
the underlying physics in a very transparent way. 

We consider a five-dimensional anti-de Sitter spacetime 
with two branes in it; a moving positive tension brane 
representing our Universe and a second brane which, for 
definitcness, is kept at rest. This setup is depicted in 
Fig. [I]. For this model we have previously shown that 
in a radiation dominated Universe, where the second, 
fixed brane is arbitrarily far away, no gravitons are 
produced pTt]. 



The particular model which we shall consider is 
strongly motivated by the ekp yrotic or cyclic Universe 
and similar ideas f| || g|, ||, H @, H§ ||. 
In this model, roughly speaking, the hot big bang 
corresponds to the collision of two branes; a moving 
bulk brane which hits "our" brane, i.e. the observable 
Universe. Within such a model, it seems to be possible to 
address all major cosmological problems (homogeneity, 
origin of density perturbations, monopole problem) 
without invoking the paradigm of inflation. For more 
details see |38 but also |59) for critical comments. 
One important difference between the ekpyrotic model 
and standard inflation is that in the latter one tensor 
perturbations have a nearly scale invariant spectrum. 
The ekpyrotic model, on the other hand, predicts a 
strongly blue gravitational wave spectrum with spectral 
tilt tit — 2 |Q. This blue spectrum is a key test 
for the ekpyrotic scenario since inflation always pre- 
dicts a slightly red spectrum for gravitational waves. 
One method to detect a background of primordial 
gravitational waves of wavelengths comparable to the 
Hubble horizon today is the polarization of the cosmic 
microwave background. Since a strongly blue spectrum 
of gravitational waves is unobservably small for large 
length scales, the detection of gravitational waves in the 
cosmic microwave background polarization would falsify 
the ekpyrotic model |3q| . 

Here we consider a simple specific model which is 
generic enough to cover important main features of the 
generation and evolution of gravitational waves in the 
background of a moving brane whose trajectory involves 
a bounce. First, the physical brane moves towards the 
static brane, initially the motion is very slow. During 
this phase our Universe is contracting, i.e. the scale 
factor on the brane decreases, the energy density on 
the brane increases and the motion becomes faster. We 
suppose that the evolution of the brane is driven by a 
radiation component on the brane, and that at some 
more or less close encounter of the two branes which we 
call the bounce, some high-energy mechanism which we 
do not want to specify in any detail, turns around the 
motion of the brane leading to an expanding Universe. 
Modeling the transition from contraction to subsequent 
expansion in any detail would require assumptions about 
unknown physics. We shall therefore ignore results 
which depend on the details of the transition. Finally 
the physical brane moves away from the static brane 
back towards the horizon with expansion first fast and 
then becoming slower as the energy density drops. This 
model is more similar to the pyrotechnic Universe of 
Kallosh, Kofman and Linde |3^| where the observable 
Universe is also represented by a positive tension brane 
rather than to the ekpyrotic model where our brane has 
negative tension. 

We address the following questions: What is the spec- 
trum and energy density of the produced gravitons, the 



massless zero mode and the KK modes? Can the gravi- 
ton production in such a brane Universe lead to limits, 
e.g. on the AdS curvature scale via the nucleosynthesis 
bound? Can the KK modes provide the dark matter 
or lead to stringent limits on these models? Similar 
results could be obtained for the free gravi-photon and 
gravi-scalar, i.e. when we neglect the perturbations of 
the brane energy momentum tensor which also couple to 
these gravity wave modes which have spin-1 respectively 
spin-0 on the brane. 

The reminder of the paper is organized as follows. 
After reviewing the basic equations of braneworld cos- 
mology and tensor perturbations in Sec. II, we discuss 
the dynamical Casimir effect approach in Sec. III. In 
Sec. IV we derive expressions for the energy density 
and the power spectrum of gravitons. Thereby we show 
that, very generically, KK gravitons cannot play the 
role of dark matter in warped braneworlds. This is 
explained by the localization of gravity on the moving 
brane which we discuss in detail. Section V is devoted to 
the presentation and discussion of our numerical results. 
In Sec. VI we reproduce some of the numerical results 
with analytical approximations and we derive fits for 
the number of produced gravitons. We discuss our main 
results and their implications for bouncing braneworlds 
in Sec. VII and conclude in Sec. VIII. Some technical 
aspects are collected in appendices. 

The main and most important results of this rather long 
and technical paper are published in the Letter |47[|. 



II. GRAVITONS IN MOVING BRANEWORLDS 

A. A moving brane in AdSs 

We consider a AdS-5 spacetime. In Poincare coordi- 
nates, the bulk metric is given by 



ds 2 = g AB dx A dx B = — \-dt 2 + Sijdx^dx 3 + dy 2 ] . 

y 

(2.1) 

The physical brane (our Universe) is located at some time 
dependent position y = yb(t), while the 2nd brane is at 
fixed position y = y s (see Fig. |l|). The induced metric on 
the physical brane is given by 



ds z 
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= a 2 (r]) [-An 2 + 5 ij dx i dx i ] , 



dt 2 + 5ijdx z dx j 



(2.2) 



is the scale factor and n denotes the conformal time of 
an observer on the brane, 



drj 



We have introduced the brane velocity 

LH 



dyb 
dt 



Vl + L 2 H 2 



and 
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7 = 



= \/l + L 2 H 2 



^/T^V 2 

Here H is the usual Hubble parameter, 

H = a/a 2 = a~ x H = -L~ V 



(2.4) 

(2.5) 
(2.6) 

(2.7) 



and an overdot denotes the derivative with respect to 
conformal time n. The bulk cosmological constant A is 
related to the curvature scale L by A = — 6/L 2 . The 
junction conditions on the brane lead to |37|, |48| 



« 5 (p + T) = 6 
k 5 ( p + P) = - 



Vl + L 2 H 2 
I ' 
2LH 

aVl + L 2 H 2 



(2.8) 
(2.9) 



Here T is the brane tension and p and P denote the 
energy density and pressure of the matter confined on 
the brane. Combining (2.8) and (2.S) results in 



-3Ha(p + P) 



(2.10) 



while taking the square of ( |2.8| ) leads to 



H 2 = -±Tp 1 
18 y{ 



2TJ 



k 2 T 2 

5 

36 



-h- r2AU 



These equations form the basis of brane cosmology and 
have been discussed at length in the literature (for re- 
views see [30|). The last equation is called the mod- 
ified Friedmann equation for brane cosmology [|l|. For 
usual matter with p + P > 0, p decreases during expan- 
sion and at sufficiently late time p -C T. The ordinary 
four-dimensional Friedmann equation is then recovered if 



<T 2 
12 



V- 



and we set 



8ttG 4 



k 2 T 
~6~ 



(2.12) 



Here we have neglected a possible four-dimensional cos- 
mological constant. The first of these equations is the 
RS fine tuning implying 



K5 = K4 L . 



(2.13) 



where 



Defining the string and Planck scales by 
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respectively, the RS fine tuning condition leads to 

2 



L 

IT 



Lpi 



(2.15) 



As outlined in the introduction, we shall be interested 
mainly in a radiation dominated low-energy phase, hence 
in the period where 

P = —p and \v\ Cl so that 7 ~ 1 ,drj ~ dt . 

In such a period, the solutions to the above equations are 
of the form 



o(t) 
Vb(t) 

«(*) 



\t\+t b 
L ' 
L 2 

|t|+*6 ' 
sgn(i)L 2 

"(l*l+^) 2 



(2.17) 
(2.18) 

(2.19) 



Negative times (t < 0) describe a contracting phase, 
while positive times (t > 0) describe radiation dominated 
expansion. At t — 0, the scale factor exhibits a kink and 
the evolution equations are singular. This is the bounce 
which we shall not model in detail, but we will have to in- 
troduce a cutoff in order to avoid ultraviolet divergencies 
in the total particle number and energy density which are 
due to this unphysical kink. We shall show, that when 
the kink is smoothed out at some length scale, the pro- 
duction of particles (KK gravitons) of masses larger than 
this scale is exponentially suppressed, as it is expected. 
The (free) parameter tb > determines the value of the 
scale factor at the bounce a b , i.e. the minimal interbrane 
distance, as well as the velocity at the bounce v b 
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a b = a(0) = — — , |v(0)| =v b = -5- 



(2.20) 



Apparently we have to demand tb > L which implies 
y b (t) < L. 



B. Tensor perturbations in AdSs 

We now consider tensor perturbations on this back- 
ground. Allowing for tensor perturbations hij(t,x,y) of 
the spatial three-dimensional geometry at fixed y, the 
bulk metric reads 



d,s 2 



T 2 

V 



+ (S i:j + 2h lj )dx l Ax J + dy 2 ] . (2.21) 



Tensor modes satisfy the traceless and transverse condi- 
tions, h\ — dihj =0. These conditions imply that hy has 
only two independent degrees of freedom, the two polar- 
ization states • = x,+. We decompose hij into spatial 



Fourier modes, 
hij(t,x,y) = 



d 3 k 
(2tt) 3 / 2 



£ e 4 *-e; i (k)/i.(t,i,;k) 



(2.22) 

where e'j (k) are unitary constant transverse-traceless po- 
larization tensors which form a basis of the two polariza- 
tion states • = x, +. For hij to be real we require 



K(t,V,^) = h.(t,y;-k). 



(2.23) 



The perturbed Einstein equations yield the equation of 
motion for the mode functions h, , which obey the Klein- 
Gordon equation for minimally coupled massless scalar 
fields in AdS, |25|, |5ll 



-0,, 



y 



h.(t,y;k)=0 



(2.24) 



In addition to the bulk equation of motion the modes also 
satisfy a boundary condition at the brane coming from 
the second junction condition, 



LHd t h. - v 7 ! + L 2 H 2 d y h. 



- 7 (vd t + d y )h.\ yb = faPIir . (2.25) 

Here n^ T) denotes possible anisotropic stress perturba- 
tions in the brane energy momentum tensor. We are in- 
terested in the quantum production of free gravitons, not 
in the coupling of gravitational waves to matter. There- 
fore we shall set II^ T> = in the sequel, i.e. we make 
the assumpti on th at the Universe is filled with a perfect 
fluid. Then, (|2.25|) reduces to 1 



(vd t + d y ) h.\ yb{t) = 



(2.26) 



This is not entirely correct for the evolution of gravity 
modes since at late times, when matter on the brane is 
no longer a perfect fluid (e.g., free-streaming neutrinos) 
and anisotropic stresses develop which slightly modify 
the evolution of gravitational waves. We neglect this sub- 
dominant effect in our treatment. (Some of the difficul- 
ties which appear when_JT T) ^ are discussed in 
The wave equa tion ( [2.24 ) together with the boundary 
condition ( 2.26 ) can also be obtained by variation of the 
action 



L 3 



S h = 2 — ^jdtjd 3 k 
x[\d t h.\ 



dy 



y b {t) V 

\dyh.\ 2 ~ k 2 \h.\ 2 



(2.27) 



which follows from the second order perturbation of the 
gravitational Lagrangian. The factor 2 in the action is 



1 In Equations (4) and (8) of our Letter |47| two sign mistakes 
have creeped in. 
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due to Z2 symmetry. Indeed, Equation ( 2.26| ) is the only 
boundary condition for the perturbation amplitude /i, 
which is compatible with the variational principle SSh = 
0, except if h, is constant on the brane. Since this issue is 
important in the following, it is discussed more detailed 
in Appendix [a[ 



Equations of motion in the late time/low energy 
limit 



In this work we restrict ourselves to relatively late 
times, when 



pT ^> (? and therefore |u| <C 1. 



(2.28) 



In this limit the conformal time on the brane agrees 
roughly with the 5D time coordinate, dr\ ~ dt and we 
shall therefore not distinguish these times; we set t = r). 
We want to study the quantum mechanical evolution 
of tensor perturbations within a canonical formulation 
similar to the dynamical Casimir effect for the electro- 
magnetic field in dynamical cavities jl5], [l6| |T^]. In or- 
der to pave the way for canonical quantization, we have 
to introduce a suitable set of functions allowing the ex- 
pansion of the perturbation amplitude h» in canonical 
variables. More precisely, we need a complete and or- 
thonormal set of eigenfunctions <p a of the spatial part 
-y 3 d y [y~ 3 d y ] of the differential opera- 



u v X u y ~ 



tor ( |2.24 ), The existence of such a set depends on the 
boundary conditions and is ensured if the problem is 
of Sturm-Liouville type (see, e.g.,[^3|). For the junc- 
tion condition ( |2.26|) , such a set does unfortunately not 
exist due to the time derivative. One way to proceed 
would be to introduce other coordinates along the lines 
of p>4| for which the junction condition reduces to a sim- 
ple Neumann boundary condition leading to a problem 
of Sturm-Liouville type. This transformation is, however, 
relatively complicated to implement without approxima- 
tions and is the subject of future work. 
Here we shall proceed otherwise, harnessing the fact that 
we are interested in low energy effects only, i.e. in small 
brane velocities. Assuming that one can neglect the 
time deriv ative in the junction condition since |v| <C 1, 
Eq. (2.25) reduces to a simple Neumann boundary con- 
dition. We shall therefore work with the boundary con- 
ditions 



d y h.\ yb = d y h.\ ys = 



(2.29) 



Then, at any time t the eigenvalu e pro blem for the spatial 
part of the differential operator ( 2.24 ) 



<t>a(t,V) 



-y 3 dy [y- 3 dy(t> a (t,y)] 

m 2 a (t)<p a (t,y) (2.30) 



is of Sturm-Liouville type if we demand that the 4>a 's are 
subject to the boundary conditions (2.29). Consequently, 



the set of eigenfunctions {<j> a (t, 2/)}^Lo is complete, 

2 Mt, v)Mt, v) = <% - y)y 3 , (2.31) 



and orthonormal with respect to the inner-product 
'' " 1 dy 



Vb 



(t) y" 



7</> a (t,y)Mt,v) = 8 a p. (2.32) 



Note the factor 2 in front of both expressions which is 
necessary in order to take the Z2 symmetry properly into 
account. 

The eigenvalues m a (t) are time-dependent and discrete 
due to the time-dependent but finite distance between 
the branes and the eigenfunctions <fi a (t,y) are time- 
dependent in particular becau se of the time dependence 
of the boundary conditions ( 2.29 ). The case a = 
with too = is the zero mode, i.e. the massless four- 
dimensional graviton. Its general solution in accordance 
with the boundary conditions is just a constant with re- 
spect to the extra dimension, <f>o{t, y) = <fro{t), and is fully 
determined by the normalization condition (0 Oj 0o) = 1 ; 



00 (*) = 



ysVbjt) 
Vy 2 s~yi(t)' 



(2.33) 



For a = i G {1, 2,3, ■■■ , | with eigenvalues m, > 0, the 
general solution of ( [2.30 ) is a combination of the Bessel 
functions J2 (Wj(t) y) and Y2 (m,i(t)y). Their particular 
combination is determined by the boundary condition at 
the moving brane. The remaining boundary condition at 
the static brane selects the possible values for the eigen- 
values TOj(t), the KK masses. For any three-momentum 
k these masses build up an entire tower of momenta in 
the j/-direction; the fifth dimension. Explicitely, the so- 
lutions (j>i(t,y) for the KK modes read 2 

<Pi{t, y) = N t (t)y 2 C 2 ( mi (t) y) (2.34) 

with 

C v (miy) = Yi(m l y b )J l ,(m i y)-Ji(m i y b )Y u (m l y). (2.35) 
The normalization reads 



Ni(t,y b ,y s 



1 



_yfCf(mi Vs) - (2/(771471-))' 
where we have used that 

2 



C 2 (77i; 2/6 ) 



7TTOi y b 



(2.36) 



(2.37) 



2 Note that we have changed the parameterization of the solutions 
with respect to |B7| for technical reasons. There, we also did not 
take into account the factor 2 related to Z2 symmetry. 
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It can be simplified further by using 

Yi(m t y b ) 2 



Ciijrii y s ) 



Yi(mi y^-nrrii y s 



leading to 



rriiTT 



Y?(m t y s ) 



Y{(miy b ) - Y{(m t y s ) 



(2.38) 



(2.39) 



Note that it is possible to have ^(mj y s ) — F^m; y b ) = 
0. But then both Y^(niiy s ) — Y^{rmyb) = and 
Eq. (2.39) has to be unde rstood as a limit. For that rea- 
son, the expression ( 2.36 ) for the normalization is used 
in the numerical simulations later on. Its denominator 
remains always finite. 

The time-dependent KK masses {mi{t)}f^ 1 are deter- 
mined by the condition 



Ci (mi(t)y s ) = 



(2.40) 



Because the zeros of the cross product of the Bessel func- 
tions Ji and Yi are not known analytically in closed form, 
t he K K-spectrum has to be determined by solving Eq. 
( |2.40 ) numerically 3 . An important quantity which we 
need below is the rate of change rhi/mi of a KK mass 
given by 



rrii 

m z = — =y b 



2 9 



N 



(2.41) 



where the rate of change of the brane motion y b is just 
the Hubble parameter on the brane 



Vb{t) 



Vb{t) 

y b {t) 



-Ha = — = —H 



(2.42) 



On account of the completeness of the eigenfunctions 
4>a(t,y) the gravitational wave am plitu de /i.(t,y;k) 
subject to the boundary conditions ( 2.29 ) can now be 
expanded as 



h.(t,y;k) 




a.k 



,(t)(f> a (t,y) 



(2.43) 



The coefficients q a ,k.*(t) are canonical variables describ- 
ing the time evolution of the perturbations and the fac- 
tor \f K5/L 3 has been introduced in order to rend er the 
Qa,k,»'s canonically normalized. In order to satisfy ( 2.23| ) 
we have to impose the same condition for the canonical 
variables, i.e. 



Qa,k,* — <7ct,-k,» 



(2.44) 



' Approximate expressions for the zeros can be found in . 



One could now insert the expansion ( 2.43| ) into the wave 
equation ( 2.24 ), multiplying it by 4>p(t, y) and integrating 
out the y— dependence by using the orthonormality to de- 
rive the equations of motion for the variables q a ,k,*- How- 
ever, as we explain in Appendix |a|, a Neumann boundary 
condition at a moving brane is not compatible with a free 
wave equation. The only consistent way to implement the 
bounda ry co nditions (2.2£) is therefore to consider the 
action ( 2.27| ) of the perturbations as the starting point 
t o de rive t he eq uations of motion for q a ,k,»- Inserting 
( [2.43 ) into ( [2.27 ) leads to the canonical action 



+ ^ [ M af3 (<?a,k,.?/3,-k,. + ?a,-k,.^,k,.) 

a/3 

+ N a/3 q a>K .qp _ k ,.] } . (2.45) 

We have introduced the time-dependent frequency of a 
graviton mode 



< fc - v/^+< , & = N , (2.46) 
and the time-dependent coupling matrices 

M a p = (d t (t>a,M , (2.47) 

N a p = {dt(t> a ,dt<t> )=^M al M 0n (2.48) 

7 

which are given explicitely in Appendix |b] (see also p7{). 
Consequently, the equations of motion for the canonical 
variables are 

+ [^ctf - N <*f] 9/3,k,« = • ( 2 - 49 ) 





The motion of the brane through the bulk, i.e. 
the expansion of the Universe, is encoded in the 
time-dependent coupling matrices M a f3,N a p. The 
mode couplings are caused by the time-dependent 
boundary condition d y h,(t,y)\ yb — which forces the 
eigenfunctions <p a (t,y) to be explicitly time-dependent. 
In addition, the frequency of a KK mode u> a> k is also 
time-dependent since the distance between the two 
branes changes when the brane is in motion. Both time- 
dependencies can lead to the amplification of tensor 
perturbations and, within a quantum theory which is 
developed in the next section, to graviton production 
from vacuum. 

Because of translation invariance with respect to the 
directions parallel to t he brane, modes with different k 
do not couple in (2.49). The three-momentum k enters 
the equation of motion for the perturbation only via the 
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frequency uj a .k, i.e. as a global quantity. Equation ( 2.49| ) 
is similar to the equation describing the time-evolution 
of electromagnetic field modes in a three-dimensional 
dynamical cavity ||l6|| and may effectively be described 
by a massive scalar field on a time-dependent interval 
|p7| . For the electromagnetic field, the dynamics of the 
cavity, or more precisely the motion of one of its walls, 
leads to photon creation from vacuum fluctuations. This 
phenomenon is usually referred to as dynamical Casimir 
effect. Inspired by this, we shall call the production of 
gravitons by the moving brane as dynamical Casimir 
effect for gravitons. 



D. Remarks and comments 

In |37j| we have already shown that in the limit where 
the fixed brane is sent off to infinity, y s — > oo, only the 
M o matrix element survives with M o = —H[l + 0(e)] 
and e = yb/y s - Moo expresses the coupling of the zero 
mode to the brane motion. Since all other couplings dis- 
appear for e — > all modes decouple from each other and, 
in addition, the canonical variables for the KK modes de- 
couple from the brane motion itself. This has led to the 
result that at late times and in the limit y s 3> y&, the KK 
modes with non- vanishing mass evolve trivially, and only 
the massless zero mode is coupled to the brane motion 
with 



9o, k,. 



H-W 



10, k, 



= 



(2.50) 



Since </> oc l/a [cf. Eqs. (|4.2|),(|4.5|)1 we have found in j37j] 
that the gravitational zero mode on the brane /io,» i~t] k) = 
y/ns/ i 3 5o,k,»0o (t, yb) evolves according to 

h ,.(t; k) + 2Hh ,.(t; k) + fc 2 h ,.(t; k) = , (2.51) 

which explicitely demonstrates that at low energies (late 
times) the homogeneous tensor perturbation equation in 
brane cosmology reduces to the four-dimensional tensor 
perturbation equation. 

An important comment is in order here concerning 
the RS II model. In the limit y s — ► oo the fixed brane 
is sent off to infinity and one ends up with a single 
positive tension brane in AdS, i.e. the RS II model. 
Even though we have shown that all couplings except 
Moo vanish in this limit, that does not imply that this 
is necessarily the case for the RS II setup. Strictly 
speaking, the above arguments are only valid in a two 
brane model with y s 3> 1. Starting with the RS II model 
from the beginning, the coupling matrices do in general 
not vanish when calculated with the corresponding 
eigenf unctions which can be found in, e.g., p2| . One 
just has to be careful when taking those limits. But 
what the above consideration demonstrates is that, if 
the couplings of the zero mode to the KK modes vanish, 
like in the y s ^> 1 limit or in the low energy RS II model 



as observed in numerical simulations (see below) the 
standard evolution equation for the zero mode emerges 
automatically from five-dimensional perturbation theory. 

Starting from five-dimensional perturbation theory, 
our formalism does imply the usual evolution equation 
for the four-dimensional graviton in a FLRW-Universe 
in the limit of vanishing couplings. This serves as a 
very strong indication (but certainly not proof!) for 
the fact that the approach based on the approximation 
( pT29[ ) and the expansion of the action in canonical 
variables rather than the wave equation is consistent 
and leads to results which should reflect the physics at 
low energies. As alr eady outlined, if one would expand 
the wave equation ( 2.24 ) in the set of functions <p a , 
the resulting equation of motion for the cor respo nding 
canonical variables is different from Eq. ( [2.49 ) and 
cannot be derived from a Lagrangian or Hamiltonian 
(see Appendix |a|). Moreover, in J|o| the low energy RS 
II scenario has been stu died numerically including the 
full junction condition (2.26) without approximations 
(see also |27]]). Those numerical results show that 
the evolution of tensor perturbations on the brane is 
four-dimensional, i.e. described by Eq. (2.51) derived 
here analytically. Combining these observations gives 
us confidence that the used approach based on the 
Neumann boundary condition approximation and the 
action as starting point for the canonical formulation is 
adequate for the study of tensor perturbations in the 
low energy limit. The many benefits this approach offers 
will become visible in the following. 



III. QUANTUM GENERATION OF TENSOR 
PERTURBATIONS 

A. Preliminary remarks 

We now introduce a treatment of quantum generation 
of tensor perturbations. This formalism is an advance- 
ment of the method which is presented in |l^, |l6|, [lTj 
for the dynamical Casimir effect for a scalar field and 
the electromagnetic field to gravitational perturbations 
in the braneworld scenario. 

The following method is very general and not restricted 
to a particular brane motion as l ong a s it complies with 
the low energy approach [cf. Eq. (2.28)]. We assume that 
asymptotically, i.e. for t — > ±oo, the physical brane 
approaches the Cauchy horizon (yb — > 0), moving very 
slowly. Then, the coupling matrices vanish an d the KK 
masses are constant (for y& close to zero, Eq. ( 2.40Q re- 
duces to Ji(rriiy s ) = 0): 



lim M a p(t) = 

t — >±oo 



lim m a (t) — const. Va,/? 

t — >±oo 



(3.1) 

In this limit, the system ( 2.49] ) reduces to an infinite set 
of uncoupled harmonic oscillators. This allows to intro- 
duce an unambiguous and meaningful particle concept, 
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i.e. notion of (massive) gravitons. 

As a matter of fact, in the numerical simulations, the 
brane motion has to be switched on and off at finite times. 
These times are denoted by t; n and t ou t , respectively. We 
introduce vacuum states with respect to times t < t m < 
and t > tout > 0. In order to avoid spurious effects influ- 
encing the particle creation, we have to chose t- m small, 
respectively t ou t large enough such that the couplings are 
effectively zero at these times. Checking the indepen- 
dence of the numerical results on the choice of t- m and 
tout guarantees that these times correspond virtually to 
the real asymptotic states of the brane configuration. 



The annihilation and creation operators {&° u £.,a° u £.} 
correspond to a meaningful definition of final state gravi- 
tons (they are associated with positive and negative fre- 
quency solutions for t > f ut) and satisfy the same com- 
mutation relations as the initial state operators. 
Initial |0, in) = |0, t < t; n ) and final |0, out) = |0, t > t ou t) 
vacuum states are uniquely defined via 4 



: J0»=0, a° U k..|t),out) = , Va, k, 



(3.9) 



The operators counting the number of particles defined 
with respect to the initial and final vacuum state, respec- 
tively, are 



B. Quantization, initial and final state 



N. 



a,k,» 



ainf sin jCrout 



a 



out f -out 
a,k,.a,k 



(3.10) 



Canonical quantization of the gravity wave amplitude 
is performed by replacing the canonical variables q a ^, 
by the corresponding operators q a ,k,» 



h,(t,y;k) = < /y| ^ ? a ,k,. (*)&*(*> v) 



(3.2) 



Adopting the Heisenberg picture to describe the quantum 
time-evol ution , it follows that <? Q ,k,» satisfies the same 
equation (2.4£) as the canonical variable q a ,k,»- 
Under the assumptions outlined above, the operator 
<?a,k,» can be written for times t < t; n as 



The number of gravitons created during the motion of 
the brane for each momentum k, quantum number a and 
polarization state • is given by the expectation value of 
the number operator N° u £ m of final-state gravitons with 
respect to the initial vacuum state |0, in): 



jyout 
J Q.k.c 



(0,in|7V°^.|0,ii 



(3.11) 



If the brane undergoes a non-trivial dynamics between 
t in < t < tout it is aa U k ..|0,in) ^ in general, i.e. 
graviton production from vacuum fluctuations takes 
place. 



1 



ft, 



(3.3) F rom Q2.22D , the expansion (§J) and Eqs.([0|), 
( |3.7|) it follows that the quantized tensor perturbation 
with respect to the initial and final state can be written 
as 



where we have introduced the initial-state frequency 

w£ fc = LU a , k (t < t in ) • (3.4) 



This expansion ensures that Eq. (2.44) is satisfied. The 
set of annihilation and creation operators {a™ k # , .} 
corresponding to the notion of gravitons for t < tm is 
subject to the usual commutation relations 



-in -inf 
a a,k,.' a a',k',.' 



-in -m 
"a,k,»i "a'.k',* 



= 5« Q '5../5 (3) (k-k') , (3.5) 



ant 



inf 



*a,k,»> a a'M\ 



0. (3.6) 



hij(t < tin,x,y) 



E 



d 3 k <t 



(2tt) 3 / 2 



2uj 



a . k 



x u' ija (t < t in , x, y, k) + h.c. 



(3.12) 



and 



hij(t > t ou t,x,y) 



E 



d 3 k a™£ 



(2tt) 3 / 2 



x u' ija (t > tout, x, y, k) + h.c. 



(3.13) 



For times t > t ou t, i-e. after the motion of the brane has 
ceased, the operator g Q ,k,» can be expanded in a similar- 
manner, 



We have introduced the basis functions 



ik x „• 



(k) 0a(*,tf). 



(3.14) 



9a,k,»(t > tout) 
1 



out 



vout t 



i u ,, t 



(3.7) which, on account of (e*j(k))* = e* (— k), satisfy 
(M* J Q (t,x,y,k))* = u' lj a {t,x,y, -k). 



with final state frequency 



J a"fc = ^a.kit > tout) 



(3.8) 



4 Note that the notations 1 , £ < ti n ) and |0,t > tout) do not mean 
that the states are time-dependent; states do not evolve in the 
Heisenberg picture. 
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Time evolution 



During the motion of the brane the time evolution of 
the field modes is des cribe d by the system of coupled 
differential equations ( |2.49[ ). To account for the inter- 
mode couplings mediated by the coupling matrix M a p 
the operator <? Qj k,» is decomposed as 



The choice of this phase for the initial condition is in 
principle arbitrary, we could as well set 9™^ = 1. But 

with this choice, e„ is independent of t ln for t < t in 
and therefore it is also at later times independent of t- la 
if only we choose t[ n sufficiently early. This is especially 
useful for the numerical work. 



<7a,k,.0) = ^2 



2w 



-in (p) / .\ | -in- 



SPY 



(*) 



0,k 



(3.15) 

The complex functions ^f\{t) also satisfy the system of 
coupl ed differential equations (2.49). With the ansatz 
( 3.15 ) the quantized tensor perturbation at any time dur- 
ing the brane motion reads 



hij(t,x, y) 

K 5 \ ~> 
L 3 2^ 
map 



(3.16) 



d 3 k af Xm 



^S(*)^>(i,x,y,k) + h.c. 



(2tt)I ^ 



2^ 



0,k 



Due to the time-dependence of the eigenfunctions <p a , 
the time-derivative of the gravity wave amplitude con- 
tains additional mode coupling contributions. Using the 
completeness and orthnormality of the <^ Q 's it is readily 
shown that 



h.(t,y;k) 



where 



Pa,-k,. (*)</><*(*> v) (3-17) 



The coupling term arises from the time dependence of 
the mode functions <p a . Accordingly, the time derivative 



hij reads 



hij(t,x,y) 



E 

.ap 



d 3 k &f tKm 

(2n)i 



2uj 



(3.19) 



p.k 



x fal(t) ^, Q (i,x,y,k)+h, 



where we have introduced the function 



(3.20) 



By comparing Eq. (3.12) and its time-derivative with 
Eqs. (3.16) and (3.19) at t = t ln one can read off the 



initial conditions for the functions e 



(/3). 



(3.21) 



(t in ) = [-iu% k 5 a (, - MpM] &% k (3.22) 



with phase 



-iwj?t, tin. 



(3.23) 



{< 



D. Bogoliubov transformations 



The two sets of annihilation and creation operators 



a,k,< 



# } and {c 



ct,k,. 



*out f 



,} corresponding to the 



notion of initial-state and final-state gravitons are re- 
lated via a Bogoliubov transformation. Ma tchin g the 
expression for the tensor perturbation Eq. ( 3.1 6| ) and 
its time-derivative Eq. ( |3.19 ) with the final state expres- 
sion Eq. (3.13) and its corresponding time-derivative at 
t — tout one finds 



X p\k 



with 



(3.24) 



•A pa, k (^out ) 



and 



U a.fe 



p,k 



"i,fc(*out) + out fallout) 



a . k 



(3.25) 



Bp a ,k{ta 



C\a\it 

^a^k 



-a.fcV L out; 



.out 

out ) 

J a,k 



(3.26) 

where we shall stick to the phase 6° u £ defined like 0„ fc 
in ( 3.23 ) for completeness. Performing the matching at 
^out = ^in the Bogoliubov transformation should become 
trivial, i.e. the Bogoliubov coefficients are subject to vac- 
uum initial conditions 



•A-aP, fe(*in) — & a p j ^aP,k{hn) — 0. 



(3.27) 



Evaluating the Bogoliubov coefficients (3.25) and (3.2£) 
for tp Ut = t m b y making use of the initial conditions 



( |3.21 ) and ( 3.22| ) shows the consi stency. Note that the 
Bogoliubov transformation (3.24) is not diagonal due 
to the inter-mode coupling. If during the motion of 
the brane the graviton field departs form its vacuum 
state one has B a p tk (t out ) 7^ 0, i.e. gravitons have been 
generated. 



By means of Eq. (|3.24| ) the number of generated 
final state gravitons (|3.11| ), which is the same for every 
polarization state, is given by 



M™ k \t>t out ) = £ <0,in|iV°£.|0,m> 
•=+,x 

= 2^|£ /3a , fe (i out )| 2 . 





(3.28) 
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Later we will sometimes interpret t out as a continuous 
variable f out — > t such that A/"" 1 ^ — * A/" a) fc(i), i- e - it 
becomes a continuous function of time. We shall call 
M a , k (t) the instantaneous particle number [see Appendix 
C2|, however, a physical interpretation should be made 
with caution. 



E. The first order system 



From the solutions of the system of differential equa- 



tions (2.49) for the complex functions e„ L the Bogoli- 
ubov coefficient B a p,k, and hence the number of cre- 
ated final state gravitons (3.28), can now be calculated. 
It is however useful to introduce auxiliary functions 
SW^SW through 



dS(*) - #£(*) + 



(/3)/ 



J a,k 



(3.29) 
(3.30) 



These are related to the Bogoliubov coefficients via 

^4/3a,fe(iout) = (3.31) 

©out* /, ,c 



/3,k 



a, re 

OUt J "I - Q. ^ out J 



out J 



(3.32) 



/ w 



A a:fc (i ut)Ci /3 fe( i out) + A+ . (^ut Wai (*out) 



where we have defined 



1 ± 



a , k 



u a ,k(t) 



(3.33) 



Using the second order differential equation for ef\ , it is 
readily shown that the functions V^kft) satisfy 



the following system of first order differential equations: 

~i [oia, fc (*^2(*)-«aa, fc (*)'/2(*)" 

E [ c a Tlfc (*)e5S(«) + <, fc (^Sw] ( 3 - 34 ) 



«2(*) 



with 



a aa,fc(*) = 



"a,fc 



1 ± 



Ja,k(t) 



J a,k 



M Q7 (t) ± -g*M 7a (t) 



7, re 



(3.36) 
(3.37) 



The vacuum initial conditions (3.27) entail the initial 
conditions 



£ {m (t- 



2S a3 e 



n m (t- 

'la,k\ t '" 



0. 



(3.38) 



With the aid of Eq. ( p. 32; ), the coefficient B a p,k(t ut), 
and therefore the number of produced gravitons, can be 
directly deduced from the solutions to this system of 
coupled first order differential equations which can be 
solved using standard numerics. 



In the next section we will show how interesting 
observables like the power spectrum and the energy den- 
sity of the amplified gravitational waves are expressed in 
terms of the number of created gravitons. The system 
(3.34, 3.35) of coupled differential equations forms the 
basis of our numerical simulations. Details of the applied 
numerics are collected in Appendix 



IV. POWER SPECTRUM, ENERGY DENSITY 
AND LOCALIZATION OF GRAVITY 

A. Perturbations on the brane 

By solving the syste m of c ouple d differential equations 
formed by Eqs. J3.34 ) and (3.35) the time evolution of 
the quantized tensor perturbation hij (t, x, y) can be com- 
pletely reconstructed at any position y in the bulk. Ac- 
cessible to observations is the imprint which the pertur- 
bations leave on the brane, i.e. in our Universe. Of par- 
ticular interest is therefore the part of the tensor pertur- 
bation which resides on the brane. It is given by eval- 
uating Eq. (2.22) at the brane position y = y b (see also 

M) 



hj(t,x,y b ) 



al 3 k 
(2tt) 3 / 2 



•el j (h)h m (t,y b ,'k) 



The motion of the brane (expansion of the Universe) en- 
ters this expres sion via the eigenfunctions (f> a (t, y b (t)). 
We shall take (4.1) as the starting point to define ob- 
servables on the brane. 

The zero- mode function <po(t) [cf. Eq. ( 2.33|)] doe s not 
depend on the extr a dim ension y. Using Eq. (2.37), one 
reads off from Eq. (2.34) that the eigenfunctions on the 
brane 4>a{t,y b ) are 



<t> a {t,vb) = y b y a {y b ) = -y a (a) (4.2) 
a 

where we have defined 




and 



Y?{m n y s ) 



Y?(m n y b ) - Y?(m n y s ) ' 



(4.3) 
(4.4) 
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for the zero- and KK modes, respectively. One immedi- 
ately is confronted with an interesting observation: the 
function y a {a) behaves differently with the expansion of 
the Universe for the zero mode a = and the KK modes 
a = n. This is evident in particular in the asymptotic 
regime y s ^> yb, i.e. yb — > (\t\,a — + oo) where, exploit- 
ing the asymptotics of Yx (see [p5|), one finds 



, ,,, , s Lirm n L lm n ir 

y (a)~l, y n {a)~ — \Yi(m n y s ) ~ -W— 

a 2 a\2y s 

(4.5) 

Ergo, 3^o is constant while J^n decays with the expansion 
of the Universe as 1/a. For large n one can approximate 
m„ ~ mr/y s and Fi(ra„y s ) ~ Fi(mr) ~ (l/ir)yj2/n [||), 
so that 



\l2na ' 



irL 2 m r , 
2y s a 2 



(4.6) 



In summary, the amplitude of the KK modes on the brane 
decreases faster with the expansion of the Universe than 
the amplitude of the zero mode. This leads to interest- 
ing consequences for the observable power spectrum and 
energy density and has a clear physical interpretation: It 
manifest the localization of usual gravity on the brane. 
As we shall show below, KK gravitons which are traces 
of the five-dimensional nature of gravity escape rapidly 
from the brane. 



B. Power spectrum 

We define the power spectrum V{k) of gravitational 
waves on the brane as in four-dimensional cosmology by 
using the rest riction of the tensor amplitude to the brane 
position (4.1): 



^P(fc)<5( 3 )(k-k') 
= (0,in /i.(t, y 6 ; k)^i(i, k') 



(4.7) 



0,in 



i.e. we consider the expectation value of the field operator 
h, with respect to the initial vacuum state at the position 
of the brane y = yb {t) ■ In order to get a physically mean- 
ingful power spectrum, averaging over several oscillations 
of the grav itat ional wave amplitude has to be performed. 
Equation (4.7) describes the observable power spectrum 
imprinted in our Universe by the four-dimensional spin-2 
graviton component of the five-dimensional tensor per- 
turbation. 

The explicit calculation of the expectation value involv- 
ing a "renormalization" of a di verg ent contribution is car- 
ried out in detail in Appendix |C 2\. The final result reads 



P(k) 



1 fc 3 



K 5 x - 



(2tt)3 L ^ 



K a , k (t)y 2 a ( 



(4.8) 



The function 1Z a ,k(t) c an be expr essed in terms of the 
Bogoliubov coefficients ( |3.25 ) and ( |3.26 ) if one considers 
tout a $ a continuous variable t: 



■>a,k\t) 



(4.9) 



Na t k{t) is t he instantaneous particle number [cf. Ap- 



pendix CI] and the function 0^ k {t) is defined in 



Eq. qcsp. 

It is important to recall that Af a ,k{t) can in general not be 
interpreted as a physical particle number. For example 
zero modes with wave numbers such that kt < 1 can- 
not be considered as particles. They have not performed 
several oscillations and their energy density cannot be 
defined in a meaningful way. 

Equivalently, expressed in terms of the complex functions 



e cfL one nn ds 



£ g(t)l 2 

in 



+ Ol k (t), (4.10) 



with O e a k given in Eq. ( |C10j ). Equation ( |4.8] ) together 



with (|1£|) or ( |4.1C| ) holds at all times. 

If one is interested in the power spectrum at early times 

kt <C 1, it is not sufficient to take only the instantaneous 



particle number Af a ,k(t) in Eq. (4.9) into account. This 



is due to the fact that even if the mode functions e^l 
are already oscillating, the coupling matrix entering the 
Bogoliubov coefficients mig ht st ill undergo a non-trivial 
time dependence [cf. Eq. ( 6.16| )1- In the next section 
we shall show explicitly, that in a radiation dominated 
bounce particle creation, especially of the zero mode, 
only stops on sub-Hubble times, kt > 1, even if the mode 
functions are plane waves right after the bounce [cf, e.g., 
Figs. ||, 0, [|. Therefore, in order to determine the per- 
turbation spectrum of the zero mod e, on e has to make 
use of the full expression expression ( 4.1C ) and may not 
use (4.11), given below. 

At late times, kt 3> 1 (t > t out ) when the brane moves 
slowly, the couplings M a p go to zero and particle cre- 
ation has come to an end, both functions 0^ k and O e a k 
do not contribute to the observable power spectrum after 
averaging over several oscillations. Furthermore, the in- 
stantaneous particle number then equals the (physically 
meaningful) number of created final state gravitons Af° n k 
and the KK masses are constant. Consequently, the ob- 
servable power spectrum at late times takes the form 



V(k,t>t out ) 



a 2 (2tt) 3 ^ u 



y 2 » , (4.ii) 



where we have used that k$/L = k±. Its dependence on 
the wave number k is completely determined by the spec- 
tral behavior of the number of created gravitons Af° u k ■ 
It is useful to decompose the power spectrum in its zero- 
mode and KK-contributions: 



V = V + Vkk- 



(4.12) 
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In the late time regime, using Eqs. (4.11) and 
zero-mode power spectrum reads 



I), the 



(4.13) 



As expected for a usual four-dimensional tensor perturba- 
tion (massless graviton) , on sub-Hubble scales the power 
spectrum decreases with the expansion of the Universe 
as 1/a 2 . 

In contrast, the KK mode power spectrum for late times, 
given by 



KK 



(k,t > t ou t) 



fc 3 KjL 2 

a 4 32tt 



(4.14) 

decreases as 1/a 4 , i.e. with a factor 1/a 2 faster than 
Vq. The gravity wave power spectrum at late times is 
therefore dominated by the zero-mode power spectrum 
and looks four dimensional. Contributions to it arising 
from five-dimensional effects are scaled away rapidly as 
the Universe expands due to the 1/a 4 behavior of "Pkk- 
In the limit of large masses m n y s > 1, n > 1 and for 
wave lengths k <C m n such that u> ni k — m n , the late-time 
KK-mode power spectrum can be approximated by 



PkkCM > *out) = 



fc 3 k 4 L 2 
a 4 16n 2 y, 



E 



KTk 



(4.15) 



where we have inserted Eq. (4J3) for y 2 (a). 
Note that the formal summations over the particle num- 
ber might be ill defined if the brane trajectory contains 
unphysical features like discontinuities in the velocity. An 
appropriate regularization is then necessary, for example, 
by introducing a physically motivated cutoff. 



C. Energy density 

For a usual four-dimensional tensor perturbation h^ v 
on a background metric an associated effective energy 
momentum tensor can be defined unambiguously by (see, 
e.g., @H) 



1 

^4 



/ill/'' 



■a/3 



(4.16) 



where the bracket stands for averaging over several pe- 
riods of the wave and denotes the covariant deriva- 
tive with respect to the unperturbed background metric. 
The energy density of gravity waves is the 00-component 
of the effective energy momentum tensor. We shall use 
the same effective energy momentum tensor to calculate 
the energy density corresponding to the four-dimensional 
spin-2 graviton component of the five-dimensional ten- 
sor perturbation on the b rane , i.e. for the perturbation 
hij(t,x,yb) given by Eq. (11). For this it is important 
to remember that in our low energy approach, and in 



particular at very late times for which we want to cal- 
culate the energy density, the conformal time rj on the 
brane is identical to the conformal bulk time t. The en- 
ergy density of four-dimensional spin-2 gravitons on the 
brane produced during the brane motion is then given by 
[see also [361 1 



1 



K4 a z 



0, in\hij(t, x, (i, x, 2,6)10, in) ) . (4.17) 



Here the outer bracket denotes averaging over several os- 
cillations, which (in contrast to the power spectrum) we 
embrace from the very beginning. The factor 1 /a 2 comes 
from the fact that an over-dot indicates the derivative 
with respect t. A detailed calculation is carried out in 
Appendix V, 31 leading to 



P=-J 



a 



d 3 k 



(4.18) 



where again Af a ,k(t) is the instantaneous particle number. 
At late times t > t ou t after particle creation has ceased, 
the energy density is therefore given by 



1 

9= -a 



E 



_fk 

(27T) 



(4.19) 



This expression looks at first sight very similar to a 
"naive" definition of energy density as integration over 
momentum space and summation over all quantum num- 
bers a of the energy u™J N° u £ of created gravitons. 
(Note that the graviton number Af° u £ already contains 



the contributions of both polarizations [see Eq. (|3.28|) 1.) 
However, the important difference is the appearance of 
the function y 2 (a) which exhibits a different dependence 
on the scale factor for the zero mode compared to the 
KK modes. 

Let us decompose the energy density into zero-mode and 
KK contributions 



P = Po + Pkk- 



(4.20) 



For the energy density of the massless zero mode one 
then obtains 



1 



Po 



d 3 fc 

(27T) 3 



(4.21) 



This is the expected behavior; the energy density of stan- 
dard four-dimensional gravitons scales like radiation. 
On contrast, the energy density of the KK modes at late 
times is found to be 



Pkk 



E 



d A k 



(2tt) 3 

(4.22) 

which decays like 1/a 6 . As the Universe expands, the en- 
ergy density of massive gravitons on the brane is there- 
fore rapidly diluted. The total energy density of gravita- 
tional waves in our Universe at late times is dominated 
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by the standard four-dimensional graviton (massless zero 
mode). In the large mass limit m n y s ^> l,n ^> 1 the KK- 
energy density can be approximated by 



which, by virtue of 



Pkk 



ttL 2 



d 3 k 
(2^)3 



(4.23) 



Due to the factor m n coming from the function y 2 , i.e. 
from the normalization of the functions <fi n (t,y), for the 
summation over the KK-tower to converge, the number 
of produced gravitons nas to decrease faster than 

1/m^ for large masses and not just faster than 1/m^ as 
one might naively expect. 



D. Escaping of massive gravitons and localization 
of gravity 

As we have shown, the power spectrum and energy 
density of the KK modes scale, at late times when par- 
ticle production has ceased, with the expansion of the 
Universe like 



Vkk oc 1/a 4 , pkk oc 1/a 6 



(4.24) 



Both quantities decay by a factor 1/a 2 faster than the 
corresponding expressions for the zero-mode graviton. In 
particular, the energy density of the KK particles on the 
brane behaves effectively like stiff matter. Mathemat- 
ically, this difference arises from the dist inct behavior 
of the functions J^o (a) and y n (a) [cf . Eq. ( |4.5| )1 and is a 
direct consequence of the warping of the fifth dimension. 
But what is the underlying physics? As we shall discuss 
now, this scaling behavior for the KK particles has 
indeed a very appealing physical interpretation which is 
in the spirit of the RS model. 

First, the mass m„ is a comoving mass. The (in- 
stantaneous) 'comoving' frequency or energy of a KK 
graviton is U! n> ). = y/k 2 + m^, with comoving wave 
number k. The physical mass of a KK mode measured 
by an observer on the brane with cosmic time dr = adt 
is therefore m n /a, i.e. the KK masses are redshiftcd 
with the expansion of the Universe. This comes from 
the fact that m n is the wave number corresponding to 
the y-direction with respect to the bulk time t which 
corresponds to conformal time r] on the brane and not to 
physical time. It implies that the energy of KK particles 
on a moving AdS brane is redshifted like that of massless 
particles. From this alone one would expect that the 
energy density of KK modes on the brane decays like 
1/a 4 (see also Appendix D of p2[). 

Now, let us define the "wave function" for a gravi- 
ton 



,.3/2 



1, satisfies 



dy*l(t,y) = 1 



(4.26) 



Fro m th e expansion of the gravity wave amplitude 
Eq. ( 2.43| ) and the normalization condition it is clear that 
^a(t, y) gives the probability to find a graviton of mass 
m a for a given (fixed) time t at position y in t he Z 2 - 
symmetric AdS-bulk. Since (j> a satisfies Equation (2.30), 
the wave function \& Q satisfies the Schrodinger like equa- 
tion 
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y 4 yl 



and the junction conditions (2.29) translate into 



(4.27) 



(4.28) 



(4.25) 



In Fig. ^| we plot the evolution of_$?(t, y) under the 
influence of the brane motion Eq. (2.18) with Vf, = 0.1. 
For this motion, the physical brane starting at y\, — ► for 
t — » — oo moves towards the static brane, corresponding 
to a contracting Universe. After a bounce, it moves 
back to the Cauchy horizon, i.e. the Universe expands. 
The second brane is placed at y s = 10L and y ranges 
from y b (t) to y s . We set ^\ = for y < y b (t) . The 
time-depe ndent KK mass mi is determined numerically 
from Eq. ( 2. 40] ). As it is evident from this Figure, ^ 2 
is effectively localized close to the static brane, i.e. the 
weight of the KK-mode wave function lies in the region 
of less warping, far from the physical brane. Thus the 
probability to find a KK mode is larger in the region 
with less warping. Since the effect of the brane motion 
on ty 2 is hardly visible in Fig. ^, we show the behavior 
of close to the physical brane in Fig. ||. This shows 
that peaks also at the physical brane but with an 
amplitude roughly ten times smaller than the amplitude 
at the static brane. While the brane, coming from 
t — ► — oo, approaches the point of closest encounter *5> 2 
slightly increases and peaks at the bounce t = where, 
as we shall show in the next Section, the production 
of KK particles takes place. Afterwards, for t — ► oo, 
when the brane is moving back towards the Cauchy 
horizon, the amplitude decreases again and so does 
the probability to find a KK particle at the position of 
the physical brane, i.e. in our Universe. The parameter 
settings used in Figures || and || are typical parameters 
which we use in the numerical simulations described 
later on. However, the effect is illustrated much better 
if the second brane is closer to the moving brane. In 
Figure k] we show for the same parameters as in 
Figures g and ^ but now with y s = L. In this case, the 
probability to find a KK particle on the physical brane 
is of the same order as in the region close to the second 
brane during times close to the bounce. However, as the 
Universe expands, rapidly decreases at the position 
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of the physical brane. 



From Eqs. (gj) and (gj) it follows that (t, y b ) cx I /a. 
The behavior of the KK-mode wave function suggests 
the following interpretation: If KK gravitons are created 
on the brane, or equivalently in our Universe, they 
escape from the brane into the bulk as the brane moves 
back to the Cauchy horizon, i.e. when the Universe 
undergoes expansion. This is the reason why the power 
spectrum and the energy density imprinted by the KK 
modes on the brane decrease faster with the expansion 
of the Universe than for the massless zero mode. 

The zero mode, on the other hand, is localized at 
the position of the moving brane. The profile of (f>o does 
not depend on the extra dimension, but the zero-mode 
wave function "Jo does. Its square is 



*o(*.v) 



vlvl 1 
/! - vl y 3 



i 



if Vs > Vb , 



such that on the brane {y = yb) it behaves as 
*o(*»2/6)-T' 



(4.29) 



(4.30) 



Equation ( 4.29| ) shows that, at any time, the zero 
mode is localized at the position of the moving brane. 
For a better illustration we show Eq. J4.29P in Fig. | 
for the same parameters as in Fig. [|. This is the 
"dynamical analog" of the localization mechanism for 
four-dimensional gravity discussed in 0. 

To establish contact with Q and to obtain a intu- 
itive phys ical d es cript ion, we rewrite the boundary value 
problem ( p7j ), ( g^ ) as a Schrodinger-like equation 



dfo a (t, y) + V{y, t)* a {y, t) = m a (t)9 a (y, t) (4.31) 



with 
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K\y\-yb{t)) 

6(\y\-y b (t)) , (4.32) 



where we have absorbed the boundary condition at the 
moving brane into the (instantaneous) volcano potential 
V(y, t) and made use of Z2 symmetry. Simi lar to the 
static case 0, at any time the potential (4.32) supports 
a sin gle bound state, the four-dimensional graviton 
(4.29), and acts as a barrier for the massive KK modes. 



The potential, ensuring localization of four-dimensional 
gravity on the brane and the repulsion of KK modes, 
moves together with the brane through the fifth dimen- 
sion. Note that with the expansion of the Universe, the 
"depth of the delta-function" becomes larger, expressing 
the fact that the localization of four-di mensi onal gravity 
becomes stronger at late times [cf. Eq. (4.30), Fig. em. 



In summary, the different scaling behavior for the 
zero- and KK modes on the brane is entirely a conse- 
quence of the geometry of the b ulk s pace-time, i.e. of 
the warping L 2 /y 2 of the metric (2J) 5 . It is simply a 
manifestation of the localization of gravity on the brane: 
as time evolves, the KK gravitons, which are traces of 
the five-dimensional nature of gravity, escape into the 
bulk and only the zero mode which corresponds to the 
usual four-dimensional graviton remains on the brane. 



This, and in particular the scaling behavior (4.24), 
remains also true if the second brane is removed, i.e. in 
the limit y s — > 00 , l eadin g to t he original RS II model. 
By looking at (4.15) and (4.23) one could at first think 
that then the KK-power spectrum and energy density 
vanish and no traces of the KK gravitons could be 
observed on the brane since both expressions behave 
as l/y s . But this is not the case since the spectrum of 
KK masses becomes continuous. In the continuum limit 
y s — > 00 the summation over the discrete spectrum to„ 
has to be replaced by an integration over continuous 
masses to in the following way: 



l£/(m„)^I f 

y s v Ti" J 



dm / (to) 



(4.33) 



/ is some function depending on the spectru m, fo r 
exam ple f (m n ) = Af° u ^- The pre-factor l/y s in ( |4.15| ) 
and (4.23) therefore ensures the existence of the proper 
continuum limit of both expressions. 
Another way of seeing this is to repeat the same 
calculations but using the eigenfunctions for the case 
with only one brane from the beginning. Those are 
5- function normalized and can be found in, e.g., [p2[ . 
They are basically the same as (2.34) except that the 
normalization is different since it depends on whether 
the fifth dimension is compact or not. In particular, on 
the brane, they have the same scale factor dependence 
as Q. 



At the end, the behavior found for the KK modes 
should not come as a surprise, since the RS II model 
has attracted lots of attention because of exactly this; 
it localizes usual four-dimensional gravity on the brane. 
As we have shown here, localization of standard four- 
dimensional gravity on a moving brane via a warped 
geometry automatically ensures that the KK modes 
escape into the bulk as the Universe expands because 
their wave function has its weight in the region of less 
warping, resulting in an KK-mode energy density on the 
brane which scales like stiff matter. 

An immediate consequence of this particular scaling 
behavior is that KK gravitons in an AdS braneworld 



Note that it does not depend on a particular type of brane motion 
and is expected to be true also in the high energy case which we 
do not consider here. 
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FIG. 2: Evolution of \??(t,y) = <f>l(t,y)/y :i corresponding to 
the probability to find the first KK graviton at time t at the 
position y in the AdS-bulk. The static brane is at y s = 10L 
and the maximal brane velocity is given by Vb = 0.1. 



FIG. 4: Evolution of *?(t,2/) for y s = L and w& = 0.1. 




FIG. 3: Evolution of as in Fig. |^ but zoomed into 

the bulk-region close to the moving brane. 



*o(*.V) 
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FIG. 5: Localization of four-dimensional gravity on a moving 



brane: Evolution of ^ (t,y) for y. 



1 and Vb = 0.1 



which should be compared with ^i(£, y) shown in Fig. [l|. 



cannot play the role of dark matter. Their energy density 
in our Universe decays much faster with the expansion 
than that of ordinary matter which is restricted to reside 
on the brane. 



V. NUMERICAL SIMULATIONS 

A. Preliminary remarks 

In this section we present results of numerical simula- 
tions fo r the bouncing model described by the equations 

In the numerical simulations we set L = 1, i.e. all 



dimensionful quantities are measured in units of the 
AdSs curvature scale. Starting at initial time t- m <C 
where the initial vacuum state |0,in) is defined, the 



system (3.34,3.35) is evolved numerically up to final time 
£ out . Thereby we set t in = ~2irN i Yl /k w ith 1 < N in £ N, 



such that 6 



o.k 



= 1 



[cf. Eq. (jOg)]. This implies 

2, i.e. independent of the three-dimensional 
momentum k a (plane wave) zero-mode solution always 
performs a fixed number of oscillations between t- m and 
the bounce at t = [cf. Eq. (3.38)]. The final graviton 



spectrum at Af^k ^ s calculated at late times t out » 1 
when the brane approaches the Cauchy horizon and 
graviton creation has ceased. This quantity is physically 
well defined and leads to the late-time power spectrum 
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(4.11) and energy density (4.19) on the brane. For 
illustrative purposes, we also plot the instantaneous 
particle number N a ,k.» (t) which also determines the 



power spectrum at all times [cf Eq.(4.9)]. In this section 
we shall use the term particle number respectively 
graviton number for both, the instantaneous particle 
number Af a ,k,*{t) as well as the final state graviton 
number Af° u £ . , keeping in mind that only the latter one 
is physically meaningful. 

There are two physical input parameters for the 
numerical simulation; the maximal brane velocity vi, (i.e. 
tb) and the position of the static brane y s . The latter 
determines the number of KK modes which fall within a 
particular mass range. On the numerical side one has to 
specify _ZVj n and tout, as well as the maximum number of 
KK modes n max which one takes into account, i.e. after 
which KK mode the system of differential equations is 
truncated. The independence of the numerical results 
on the choice of the time parameters is checked and the 
convergence of the particle spectrum with increasing 
"■max is investigated. More detailed information on 
numerical issues including accuracy considerations are 
collected in Appendix 



One strong feature of the brane motion (2.18) is its kink 
at the bounce t = 0. In order to study how particle 
production d epend s on the kink, we shall compare 
the motion ( 2.18 ) with the following motion which 
has a smooth transition from contraction to expansion 

(L = iy. 



Vbit) 



{\t\ + t b -t s )- 1 

a + {b/2)t 2 + (c/4)i 4 



|i| > t s 
|i| < t s 



(5.1) 



with the new parameter t s in the range < t s < This 
motion is constructed such that its velocity at \t\ = t s 
is the same as the velocity of the kink motion at the 
bounce. This will be the important quantity determin- 
ing the number of produced gravitons. For t s — > the 
motion with smooth transition approaches (2.18). The 
parameters a, b and c are obtained by matching the mo- 
tions and the first and second derivatives. Matching also 
the second derivative guarantees that possible spurious 
effects contributing to particle production are avoided. 
The parameter t s has to be chosen small enough, t s <C 1, 
such that the maximal velocity of the smooth motion is 
not much larger than Vb in order to have comparable sit- 
uations. 

For reasons which will become obvious in the next two 
sections we shall discuss the cases of long k <C 1 and 
short wavelengths k >■ 1, separately. 



B. Generic results and observations for long 
wavelengths fc< 1 

Figure ^ displays the results of a numerical simula- 
tion for three-momentum k = 0.01, static brane position 
y s = 10 and maximal brane velocity v b — 0.1. Depicted is 
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time t 

FIG. 6: Evolution of the graviton number Af a ,k,m(t) for the 
zero mode and the first ten KK modes for three-momentum 
k = 0.01 and v b = 0.1, y s = 10. 
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FIG. 7: Nn,k,»{i) for the zero mode and the first ten KK 
modes for the parameters of Fig. ^j, but without coupling of 
the zero mode to the KK modes, i.e. Mio = 0. 



the graviton number for one polarization Af a ,k,»(t) for the 
zero mode and the first ten KK modes as well as the evo- 
lution of the scale factor a(t) and the position of the phys- 
ical brane yb(t). Initial and final times are 7Vj n = 5 and 
t out = 2000, respectively. The KK-particle spectrum will 
be discussed in detail below. One observes that the zero- 
mode particle number increases slightly with the expan- 
sion of the Universe towards the bounce at t = 0. Close 
to the bounce A/o,fe,«(t) increases drastically, shows a lo- 
cal peak at the bounce and, after a short decrease, grows 
again until the mode is sub- horizon (kt 3> 1). Inside the 
horizon Ao,fe,»(i) is oscillating around a mean value with 
diminishing amplitude. This mean value which is reached 
asymptotically for t — > oo corresponds to the number of 
generated final state zero-mode gravitons A/q . Produc- 
tion of KK-mode gravitons takes effectively place only at 
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the bounce in a step-like manner and the graviton num- 
ber remains constant right after the bounce. 
In Fig. ^ we show the numerical results obtained for the 
same parameters as in Fig. [|but without coupling of the 
zero mode to the KK modes, i.e. Mio = (and thus 
also Nio = Noi = 0). One observes that the production 
of zero-mode gravitons is virtually not affected by the 
artificial decoupling 6 . Note that even if M$j = (see 
Eqs. B2), which is in general true for Neumann bound- 
ary conditions, the zero mode <?o.k.» couples in Eq. ( 2.49| ) 



four-dimensional graviton for long wavelengths is not in- 
fluenced by the KK modes. Thereby the power spectrum 
is shown before and after ave raging over several oscilla- 
tions, i.e. employing Eq. (L9) with and without the term 



to the KK modes via N j = M 00 Mj and through the 
anti-symmetric combination M a p ~ Mp a . 
In contrast, the production of the first ten KK modes is 
heavily suppressed if M,o = 0. The corresponding final- 
state graviton numbers Af° n ^ . are reduced by four orders 
of magnitude. This shows that the coupling to the zero 
mode is essential for the production of massive gravitons. 
Later we will see that this is true for light KK gravitons 
only. If the KK masses exceed m, ~ 1, they evolve in- 
dependently of the four-dimensional graviton and their 
evolution is entirely driven by the intermode couplings 
Mij. It will also turn out that the time-dependence of 
the KK mass rrii plays only an inferior role for the gen- 
eration of massive KK modes. On the other hand, the 
effective decoupling of the evolution of the zero mode 
from the KK modes occurs in general as long as k -C 1 
is satisfied, i.e. for long- wavelengths. We will see that it 
is no longer true for short wavelengths k 1. 
The effective decoupling of the zero-mode evolution from 
the KK modes makes it possible to derive analytical ex- 
pressions for the number of zero-mode gravitons, their 
power spectrum and energy density. The calculations 
are carried out in section VI A| 

In summary we emphasize the important observation 
that for long wavelengths the amplification of the four 
dimensional gravity wave amplitude during the bounce 
is not affected by the evolution of the KK gravitons. We 
can therefore study the zero mode separately from the 
KK modes in this case. 



C. Zero mode: long wavelengths k <C 1 

In Figure || we show the numerical results for the num- 
ber of generated zero-mode gravitons A/o,fc,.(t) and the 
evolution of the corresponding power spectrum Vo(k) on 
the brane for momentum k = 0.01, position of the static 
brane y s = 10 and maximal brane velocity v b = 0.1. The 
results have been obtained by solving the equations for 
the zero mode alone, i.e. without the couplings to the KK 
modes, since, as we have just shown, the evolution of the 



G{f k , respectively. Right after the bounce where the gen 
eration of gravitons is initiated and which is responsible 
for the peak in Afa,k,» at t — 0, the number of gravitons 
first decreases again. Afterwards Afo,k,* grows further un- 
til the mode enters the horizon at kt = 1. Once on sub- 
horizon scales kt ^> 1, the number of produced gravitons 
oscillates with a diminishing amplitude and asymptoti- 
cally approaches the final state graviton number ■A/g'jf.. 
During the growth of J\fo,k,» after the bounce, the power 
spectrum remains practically constant. Within the range 
of validi ty it i s in good agreement with the analytical pre- 
diction fl6.22| ) yielding (L 2 (2n) 3 / K 4 )V (k , t) = Av b {kL) 2 . 
When pa rticl e creation has ceased, the full power spec- 
trum Eq.(4.8) starts to oscillate with an decreasing am- 
plitude. The time-averaged power spectrum obtained 
by using Eq. ( |4.9| ) without the 0^,-term is perfectly 
in agreement with the analytical expression Eq. ( 6.20| ) 
which gives (L 2 (2tt) 3 / K A )V {k,t) = 2v b /t 2 . Note that at 
early times, the time-averaged power spectrum behaves 
not in the same way as the full one, demonstrating the 
importance of the term 0^ k . 

Figure |^ shows a summary of numerical results for the 
number of created zero- mode gravitons J\fo,k,»(t) for dif- 
ferent values of the three-momentum k. The maximum 
velocity at the bounce is v b — 0.1 and the second brane is 
at y s = 10. These values are representative. Other values 
in accordance with the considered low-energy regime do 
not lead to a qualitatively different behavior. Note that 
the evolution of the zero mode does virtually not depend 
on the value of y s as long as y s » 2/b(0) (see below). Ini- 
tial and final integration times are given by iVj n = 5 and 
^out = 20000, respectively. 

For sub-horizon modes we compare the fi nal graviton 
spectra with the analytical prediction (6.17). Both are in 
perfect agreement. On super-horizon scales where parti- 
cle creation has not ceased yet Afo.k,* is independent of 
k. The corresponding time-evolution of the power spec- 
tra Vo{k, t) is depicted in Fig. [l(]. For the sake of clarity, 
only the results for t > 0, i.e. after the bounce, are shown 
in both figures. 

The numerical simulations and the calculations of sec- 
tion VIA reveal that the power spectrum for the four- 



dimensional graviton for long wavelengths is blue on 
super-horizon scales, as expected for an ekpyrotic sce- 
nario. 



Quantitatively it is M Q . k ..{t = 2000) = 965.01 with and 
A/o,fc,.(* = 2000) = 965.06 without M l0 . Note that this differ- 
ence lies indeed within the accuracy of our numerical simulations 
(see Appendix |l5| ) 



The analytical calculations performed in section VIA 
rely on the assumption that y b <C y s and ti n — > — oo. 
Figure |ll| shows the behavior of the number of generated 
zero-mode gravitons of momentum k = 0.01 in depen- 
dence on the inter-brane distance and the initial integra- 
tion time. The brane velocity at the bounce is v b = 0.1 
which implies that at the bounce the moving brane is 
at yb(0) — \/vb — 0.316 (L = 1). In case of a close en- 
counter of the two branes as for y s — 0.35, the production 
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FIG. 8: Time evolution of the number of created zero-mode 
grav itons J\fo,k,»(t) and of the zero-mode power spectrum 
(Ejj): (a) for the entire integration time; (b) for t > only. 
Parameters are k — 0.01, y s = 10 and Vb = 0.1. Initial and fi- 
nal time of integration are given by iVi n = 10 and t out = 4000, 
respectively. The power spectrum is shown with and without 
the term Ooffc,., i- e - before and after averaging, respectively, 
and compared with the analytical results. 
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FIG. 9: Numerical results for the time evolution of the num- 
ber of created zero- mode gravitons Ao,fc,.(t) after the bounce 
t > for different three-momenta k. The maximal brane 
velocity at the bounce is Vb = 0.1 and the second brane is 
positioned at y 3 = 10. In the final particle spectrum the nu- 
meri cal v alues are compared with the analytical prediction 
Eq. ( |S,17| ), Initial and final time of integration are given by 
Mn = 5 and t out = 20000, respectively. 



FIG. 10: Evolution of the zero-mode power spectrum after the 
bounce t > corresponding to the values and parameters of 
Fig. ^. The num erical result s are compared to the analytical 
predictions Eqs. (6.2C) and (3.22). 
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FIG. 11: Dependence of the zero-mode particle number on 
inter-brane distance and initial integration time for momen- 
tum k = 0.01, maximal brane veloc ity V h = 0.1 in comparison 
with the analytical expression Eq. ( p,16|) . (a) Evolution of the 
instantaneous particle number TVo.fc,. (t) with initial integra- 
tion time given by iVj n = 5 for y s = 0.35, 0.5 and 1. (b) Final 
zero-mode graviton spectrum Afo,k,*(t ut = 2000) for various 
values of y s and Ni B . (c) Close-up view of (b) for large y s . 



of massless gravitons is strongly enhanced compared to 
the analytical result. But as soon as y s > 1, (i.e. y s > L) 
the numerical result is ve ry well described by the analyt- 
ical expression Eq. ( 6.1 6| ) derived under the assumption 
Vs 3> yb- For y s > 10 the agreement between both is 
very good. From panels (b) and (c) one infers that the 
numerical result becomes indeed independent of the ini- 



tial integration time when increasing N- lu . Note that in 
the limit A^i n 3> 1 the numerical result is slightly larger 
than the analytical prediction but the difference between 
both is negligibly small. This confirms the correctness 
and accur acy of the analytical expressions derived in Sec- 
tion 



VI A for the evolution of the zero-mode graviton. 
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FIG. 12: Final state KK-graviton spectra for k — 0.001, y s = 
100, different maximal brane velocities Vb and N[ n = 1, tout = 
400. The nume rical results are compared with the analytical 
prediction Eq. (5.34) (dashed line). 



FIG. 13: Final state KK-graviton spectra for k — 0.01, y 3 = 
100, different Vb and Ni n = 1, i ou t = 400. The numerical 



results are compared with the analytical prediction Eq. (6.34) 
(dashed line). For Vb = 0.3,0.5 the spectra obtained without 
KK-intermode and self-couplings (My = Vi,j) are shown 
as well. 



D. Kaluza-Klein-modes: long wavelengths k <C 1 



Because the creation of KK gravitons ceases right 
after the bounce [cf Fig. || one can stop the numerical 
simulation and read out the number of produced KK 
gravitons Af° u u . at times for which the zero mode is still 
super-horizon. 

Even though Eq. (2.4C) cannot be solved analytically, 
the KK masses can be approximated by m„ ~ mr/y s . 
This expression is the better the larger the mass. 
Consequently, for the massive modes the position of the 
second brane y s determines how many KK modes belong 
to a particular mass range Am. 

In Figure [l2| we show the KK-graviton spectra J\f° u £ . for 
three-momentum k — 0.001 and second brane position 
y s = 100 for maximal brane velocities v b = 0.1,0.3 
and 0.5. For any velocity Vb two spectra obtained 



with n r 



60 and 80 KK modes taken into account 



in the simulation are compared to each other. This 
reveals that the numerical results are stable up to a 
KK mass m„ ~ 1. One infers that first, Af° u £ . grows 
with increasing mass until a maximum is reached. The 
position of the maximum shifts slightly towards larger 
masses with increasing brane velocity vj,. Afterwards, 
Af° u u . declines with growing mass. Until the maximum 
is reached, the numerical results for the KK-particle 
spectrum are very stable. This already indicates that 
the KK-intcrmode couplings mediated by Mm are not 
very strong in this mass range. In Figure O we show 
the final KK-particle spectrum for the same parameters 
as in Fig. 03 but for three-momentum k — 0.01 and 



the additional velocity Vb — 0.9 7 . We observe the same 
qualitative behavior as in Fig. In addition we show 
numerical results obtained for v b = 0.3 and 0.5 without 
the KK-intermode and self couplings, i.e. we have set 
M^ = Vi,j by hand. One infers that for KK masses, 
depending slightly on the velocity Vb but at least up to 
m n ~ 1, the numerical results for the spectra do not 
change when the KK-intermode coupling is switched off. 
Consequently, the evolution of light, i.e. m n < 1, KK 
gravitons is virtually not affected by the KK-intermode 
coupling. 

In addition we find that also the time-dependence of the 
KK masses is not important for the production of light 
KK gravitons which is explicitly demonstrated below. 
Thus, production of light KK gravitons is driven by 
the zero-mode evolution only. Th is allows us to find 
an analytical expression, Eq. ( 6.34| ), for the number of 
produced light KK gravitons in terms of exponential 
integrals. The calculations which are based on several 



approximations are performed in Section VI C . 
In Figs. [IJ and |l3| the analytical prediction ( 6.34 ) for 



the spectrum of final state gravitons has already been 
included (dashed lines). Within its range of validity 
it is in excellent agreement with the numerical results 
obtained by including the full KK-intermode coupling. 
It perfectly describes the dependence of Af° u ^ . on the 
three-momentum k and the maximal velocity Vb- For 
small velocities Vb < 0.1 it is also able to reproduce the 
position of the maximum. This reveals that the KK- 



7 Such a high brane velocity is of c ourse not consistent with a 
Neumann boundary condition Eq. (2.29) at the position of the 
moving brane. 
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intermode coupling is negligible for light KK gravitons 
and that their production is entirely driven by their 
coupling to the four-dimensional graviton. 
The analytical prediction is very precious for testing the 
goodness of the parameters used in the simulations, in 
particular the initial time t m (respectively Ni n ). Since it 
has been derived for real asymptotic initial conditions, 
tin - * — oo, its perfect agreement with the numerical 
results demonstrates that the values for N- m used in the 
numerical simulations are large enough. No spurious 
initial effects contaminate the numerical results. 
Note, that the numerical values for J\f° u £ t in the ex- 
amples shown are all smaller than one. However, for 
smaller values of k than the ones which we consider here 
for purely numerical reasons, the number of generated 
KK-mode particles is enhanced since Af° u ^ . oc 1/fc as 
can be inferred from Eq. ( 6.34 ) in the limit k <C m n . 



If we go to smaller values of y s , fewer KK modes 
belong to a particular mass range. Hence, with the same 
or similar number of KK modes as taken into account 
in the simulations so far, we can study the behavior 
of the final particle spectrum for larger masses. These 
simulations shall reveal the asymptotical behavior of 
■Wnk» f° r m ™ ~^ 00 ano - therefore the behavior of the 
total graviton number and energy density. Due to the 
kink in the brane motion we cannot expect that the 
energy density of produced KK-mode gravitons is finite 
when summing over arbitrarily high frequency modes. 
Eventually, we will have to introduce a cutoff setti ng the 
sc ale a t which the kink-approximation [cf. Eqs. (2.17) 
- ( 2.19 )1 is 110 longer valid. This is the scale where the 
effects of the underlying unspecified high-energy physics 
which drive the transition from contraction to expansion 
become important. The dependence of the final particle 
spectrum on the kink will be studied later on in this 
section in detail. 
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FIG. 14: Final state KK-graviton spectra for k — 0.01, y 3 = 
10, different maximal brane velocities Vb and Ni n = 2, t out = 
400. The nume rical results are compared with the analytical 
prediction Eq. (|6.34[) (dashed line). 
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In Figures [14] and [15| we show final KK-graviton 
spectra for y s = 10 and three-momentum k = 0.01 and 
k = 0.1. The analytical expression Eq. (6.34) is depicted 



as well and the spectra are always shown for at least 
two values of ra max to indicate up to which KK mass 
stability of the the numerical results is guaranteed. 
Now, only two KK modes are lighter than m = 1. For 
these modes the analytical expression Eq. J6.34 ) is valid 
and in excellent agreement with the numerical results, in 
particular for small brane velocities vi, ~ 0.1. As before, 
the larger the velocity Vb the more visible is the effect of 
the truncation of the system of differential equations at 

For k — 0.01 the spectrum seems to follow a power law 
decrease right after the maximum in the spectra. In 
case of Vb = 0.1 the spectrum is numerically stable up to 
masses m n ~ 20. In the region 5 < m n < 20 the spec- 
trum is very well fitted by a power law AT° U £ . oc m~ 2 ' 7 . 
Also for larger velocities the decline of the spectrum 
is given by the same power within the mass ranges 



FIG. 15: Final state KK-graviton spectra for k — 0.1, y a — 10, 
different maximal brane velocities Vb and N- m = 2, tout = 
400. The nume rical results are compared with the analytical 
prediction Eq. (|6.34|) (dashed line). 



where the spectrum is numerically stable. For k = 0.1, 
however, the decreasing spectrum bends over at a mass 
around m„ ~ 10 towards a less steep decline. This is 
in particular visible in the two cases with Vb — 0.1 and 
0.3 where the first 100 KK modes have been taken into 
account in the simulation. The behavior of the KK-mode 
particle spectrum can therefore not be described by a 
single power law decline for masses m n > 1. It shows 
more complicated features instead, which depend on the 
parameters. We shall demonstrate that this bending 
over of the decline is related to the coupling properties of 
the KK modes and to the kink in the brane motion. But 
before we come to a detailed discussion of these issues, 
let us briefly confront numerical results of different y s to 
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FIG. 16: Upper panel: Final state KK-particle spectra for 
k = 0.01, v b = 0.1 and diff erent y s = 3,10,30 and 100. 
The analytical prediction Eq. (6.34) is shown as well (dashed 
line). Lower panel: Energy M^kAfnk • °f the produced fi- 
nal state gravitons binned in mass intervals Am = 1 for 
y s = 10,30,100. 
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0.1. 



Upper panel: Final state KK-particle spectra for 
Vb — 0.1 and diff erent y s — 3, 10, 30 and 100. 



The analytical prediction Eq. (6.34) is shown as well (dashed 
line). Lower panel: Energy uj^^Af° n £ . of the produced fi- 
nal state gravitons binned in mass intervals Am = 1 for 
y s = 10, 30, 100. 



demonstrate a scaling behavior. 

In the upper panel of Figures [l6] and ^ we com- 
pare the final KK-spectra for several positions of the 
second brane y s — 3, 10, 30 and 100 obtained for a 
maximal brane velocity Vb — 0.1 for k — 0.01 and 
0.1, respectively. One observes that the shapes of the 
spectra are identical. The bending over in the decline 
of the spectrum at masses m„ ~ 1 is very well visible 
for k = 0.1 and y s = 3, 10. For a given KK mode n 
the number of particles produced in this mode is the 
larger the smaller y s . But the smaller y s , the less KK 
modes belong to a given mass interval Am. The energy 
transferred into the system by the moving brane, which 
is determined by the maximum brane velocity Vb, is the 
same in all cases. Therefore, the total energy of the pro- 
duced final state KK gravitons of a given mass interval 
Am should also be the same, independent of how many 
KK modes are contributing to it. This is demonstrated 
in the lower panels of Figs. [l6] and |l~7| where the energy 
w nt^nl'«(' n units of L) of the generated KK gravitons 
binned in mass intervals Am = 1 is shown 8 . One 
observes that, as expected, the energy transferred into 
the production of KK gravitons of a particular mass 
range is the same (within the region where the numerical 
results are stable), independent of the number of KK 
modes lying in the interval. This is in particular evident 
for y s = 30, 100. The discrepancy for y s = 10 is due 
to the binning. As we shall discuss below in detail, the 



particle spectrum can be split into two different parts. 
The first part is dominated by the coupling of the zero 
mode to the KK modes (as shown above), whereas the 
second part is dominated by the KK-intermode couplings 
and is virtually independent of the wave number k. 
As long as the coupling of the zero mode to the KK 

modes is the dominant contribution to KK-particle 

production it is Af™t, oc l/k [cf. Eq. fl6.34|)]. Hence, 



cout 
°n.k,i 



The energy for the case y s = 3 is not shown because no KK mode 
belongs to the first mass interval. 



oc l/k if m„ > k. This explains why 
the energy per mass interval Am is one order larger for 
k = 0.01 (cf Fig. ffi| than for k = 0.1 (cf Fig. [j]]) . 

Let us now discuss the KK-spectrum for large masses. 
The qualitative behavior of the spectrum Af° n £ . and the 
mass at which the decline of the spectrum changes are 
independent of y s . This is demonstrated in Figure |l8| 
where KK-spectra for Vb — 0.1, k = 0.1, y s = 10 [cf 
Fig. [ll| and y s = 3 [cf Fig. |l7j are shown. The results 
obtained by taking the full intermode coupling into 
account are compared to results of simulations where 
we have switched off the coupling of the KK modes to 
each other as well as their self-coupling (My = Mi, j). 
Furthermore we display the results for the KK-spectrum 
obtained by taking only the KK-intermode couplings 
into account, i.e. M^o = M« = Vi. One infers that 
for the lowest masses the spectra obtained with all 
couplings are identical to the ones obtained without 
the KK-intermode (My = 0, i ^ j) and self-couplings 
(Ma = 0). Hence, as already seen before, the primary 
source for the production of light KK gravitons is their 
coupling to the evolution of the four-dimensional gravi- 
ton. In this mass range, the contribution to the particle 
creation coming from the KK-intermode couplings is 
very much suppressed and negligibly small. 
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FIG. 18: KK-particle spectra for three-momentum k = 0.1, 
maximum brane velocity vb = 0.1 and y s = 3 and 10 with 
different couplings taken into account. Th e das hed lines indi- 
cates again the analytical expression Eq. (|3.34|). 



FIG. 19: Comparison of KK-particle spectra for y 3 = 3, 
Vb = 0.1 and three-momentum k = 0.01, 0.03, 0.1 and 1 
demonstrating the independence of the spectrum on k for 
large masses. n max = 60 KK modes have been taken into 
account in the simulations. 



For masses m n ~ 4 a change in the decline of the 
spectrum sets in and the spectrum obtained without 
the coupling of the KK modes to the zero mode starts 
to diverge from the spectrum computed by taking all 
the couplings into account. While the spectrum without 
the KK-intermode couplings decreases roughly like a 
power law J\f° u £ # oc m~ 3 the spectrum corresponding 
to the full coupling case changes its slope towards 
a power law decline with less power. At this point 
the KK-intermode couplings gain importance and the 
coupling of the KK modes to the zero mode looses 
influence. For a particular mass m c — 9 the spectrum 
obtained including the KK-intermode couplings only, 
crosses the spectrum calculated by taking into account 
exclusively the coupling of the KK modes to the zero 
mode. After the crossing, the spectrum obtained by 
using only the KK-intermode couplings approaches the 
spectrum of the full coupling case. Both agree for large 
masses. Thus for large masses m n > m c the production 
of KK gravitons is dominated by the couplings of the 
KK modes to each other and is not influenced anymore 
by the evolution of the four-dimensional graviton. This 
crossing defines the transition between the two regimes 
mentioned before: for masses m„ < m c the production 
of KK gravitons takes place due to their coupling to 
the zero mode M,o, while it is entirely caused by the 
intermode couplings M^- for masses m n > m c . 
Decoupling of the evolution of the KK modes from 
the dynamics of the four-dimensional graviton for 
large masses implies that KK-spectra obtained for the 
same maximal velocity are independent of the three- 
momentum k. This is demonstrated in Fig. |l9| where we 
compare spectra obtained for = 0.1 and y s = 3 but 
different k. As expected, all spectra converge towards 
the same behavior for masses m n > m c . 
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FIG. 20: KK-particle spectra for three-momentum k = 0.1, 
maximum brane velocities Vb = 0.1 and y s = 3 for n max = 40 
obtained for different coupling combinations. 



Figure. |2C| shows KK-particle spectra for k = 0.1, = 0.1 
and y s = 3 obtained for different couplings. This plot 
visualizes how each particular coupling combination 
contributes to the production of KK gravitons. It 
shows, as already mentioned before but not shown 
explicitly, that the Ma coupling which is the r ate o f 
change of the corresponding KK mass [cf. Eqs. ( 2.41 ) 
and (B4)] is not important for the production of KK 
gravitons. Switching it off does not affect the final 
graviton spectrum. We also show the result obtained 

4- t . \ 



with all couplings but with a% (t) = u> - n fc and a u (t) = 



i.e. the time-dependence of the frequency [cf. Eq. ( [3.36 )] 
has been neglected. One observes that in this case the 
spectrum for larger masses is quantitatively slightly 
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FIG. 21: KK-particle spectra for y s = 10, Vb = 0.1, n max = 
100 and three-momentum k = 0.01 and 0.1 with different 
coup lings taken into account. The thin dashed lines indicates 
Eq. fj6 34j) and the thick dashed line Eq. (E3). 



FIG. 22: KK-particle spectra for three-momentum fc = 
0.1, i/ s = 10 and maximum brane velocities Vb = 0.03,0.1 and 
0.3 with n max = 100. As in Fig. Ell different couplings 
have been taken into account and thin aashed lines indicates 
Eq. (]6.34|) and the thick dashed line Eq. (O). 



different but has a identical qualitative behavior. If, 
on the other hand, all the couplings are switched off 
M a p = Va,(9 and only the time-dependence of the 
frequency u^fc is taken into account, the spectrum 
changes drastically. Not only the number of produced 
gravitons is now orders of magnitude smaller but also 
the spectral tilt changes. For large masses it behaves as 



Consequently, the time-dependence of 



the graviton frequency itself plays only an inferior role 
for production of KK gravitons. 

The bottom line is that the main sources of the produc- 
tion of KK gravitons is their coupling to the evolution of 
the four-dimensional graviton (May) and their couplings 
to each other (My, i ^ j) for small and large masses, 
respectively. Both are caused by the time-dependent 
boundary condition. The time-dependence of the oscilla- 
tor frequency cjj^ — mj(t) + k 2 is virtually irrelevant. 

Note that this situation is very different from ordinary 
inflation where there are no boundaries and particle 
production is due entirely to the time dependence of the 
frequency 9 . 

The behavior of the KK-spectrum, in particular 
the mass m c at which the KK-intermode couplings start 
to dominate over the coupling of the KK modes to 
the zero mode depends only on the three-momentum 
k = |k| and the maximal brane velocity Vh- This is now 
discussed. In Figure 21 we show KK-particle spectra 
for y s = 10, Vf, — 0.1, rt max = 100 and three-momenta 
k = 0.01 and 0.1. Again, the spectra obtained by taking 



Note, however, that the time-dependent KK mass mj(t) enters 
the intermode couplings. 



all the couplings into account are compared to the case 
where only the coupling to the zero mode is switched 
on. One observes that for k = 0.01 the spectrum is 
dominated by the coupling of the KK modes to the 
zero mode up to larger masses than it is the case for 
k = 0.1. For k = 0.01 the spectrum obtained taking 
into account Afjo only is identical to the spectrum 
obtained with the full coupling up to m n ~ 10. In case 
of k = 0.1 instead, the spectrum is purely zero mode 
dominated only up to m n ~ 5. Hence, the smaller the 
three-momentum k the larger is the mass range for 
which the KK-intermode coupling is suppressed, and 
the coupling of the zero mode to the KK modes is the 
dominant source for the production of KK gravitons. As 
long as the coupling to the zero mode is the primary 
source of particle production, the spectrum declines 
with a power law oc m^ 3 . Therefore, in the limiting case 
k — * when the coupling of the zero mode to the KK 
modes dominates particle production also for very large 
masses it is ■A/"°g lifc _ >0) . oc l/m 3 . 

Figure ^2| shows KK-graviton spectra obtained for 
the same parameters as in Fig. ||l] but for fixed k = 0.1 
and different maximal brane velocities Vh. Again, the 
spectra obtained by taking all the couplings into account 
are compared with the spectra to which only the coupling 
of the KK modes to the zero mode contributes. The 
mass up to which the spectra obtained with different 
couplings are identical changes only slightly with the 
maximal brane velocity Vb- Therefore, the dependence 
of m c on the velocity is rather weak even if Vb is changed 
by an order of magnitude, but nevertheless evident. 
This behavior of the spec trum can indeed be understood 
qualitatively. In Section VIC we demonstrate that the 



coupling strength of the KK modes to the zero mode 
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FIG. 23: KK-particle spectra for three-momentum k = 0.01, 
0.03, 0.1 and 1 for y s = 3 and maximum brane velocity 
Vb = 0.1 with different coup ling s taken into account where 
the notation is like in Fig. E2. From the crossing of the 
Ma — Mij — 0- and Mu = Mm = results we determine 
the fc-dep ende nce of m c (k,vt). The thick dashed line indi- 
cates Eq. (5.4). 



at the bounce t = 0, where production of KK gravitons 
takes place, is proportional to 



(5.2) 



Vb 

k 



The larger this term the stronger is the coupling of 
the KK modes to the zero mode, and thus the larger 
is the mass up to which this coupling dominates over 
the KK-intermode couplings. Consequently, the mass 
at which the tilt of the KK-particle spectrum changes 
depends strongly on the three-momentum k but only 
weakly on the maximal brane velocity due to the square 
root behavior of the coupling strength. This explains 
qualitatively the behavior obtained from the numerical 
simulations. 

An approximate expression for m c (k,Vb) can be 
obtained from the numerical simulations. In Figure |23| 
we depict the KK-particle spectra for three-momentum 
k = 0.01, 0.03, 0.1 and 1 for y s = 3 and maximum brane 
velocity Vb = 0.1 with different couplings taken into 
account. The legend is as in Fig. |2^. From the crossings 
of the — 0, i ^ j and Mu — Mio = results one 
can determine the /c-dependence of m c . Note that the 
spectra are not numerically stable for large masses, but 
they are stable in the range where m c lies [cf., e.g., 
Fig. H, for k = 0.1]. Using the data for k = 0.01,0.03 
and 0.1 one finds m c (k,Vb) oc 1/Vk . 
In Fig. ^4] KK-graviton spectra are displayed for 
k = 0.1, y s = 3 and maximal brane velocities 
v h = 0.3,0.2,0.1,0.08,0.05 and 0.03 with different 
couplings taken into account. It is in principle possible 
to determine the i^-dependence of m c from the crossings 
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FIG. 24: KK-graviton spectra for three-momentum k = 
0.1, y s = 3 and maximum brane velocities Vb = 
0.3, 0.2, 0.1, 0.08, 0.05 and 0.03 with different couplings taken 
into account where the notation is like in Fig. E2J. From the 
crossing of the Mu = Mij = 0- and Ma = Mm = results we 
determine the ^-dependence of m c . 



of the Mij =0, i 7^ j- and Mu = Mio = results as 
done for the ^-dependence. However, the values for m c 
displayed in the Figures indicate that the dependence 
of m c on Vb is very weak. From the given data it is not 
possible to obtain a good fitting formula (as a simple 
power law) for the ^-dependence of m c . (In the range 
0.1 < Vb < 0.3 a very good fit is m c = 1.12-kv® 13 / Vk.) 
The reason is twofold. First of all, given the complicated 
coupling structure, it is a priori not clear that a simple 
power law dependence exists. Recall that also the ana- 



lytical expression for the particle number Eq. (3.34) has 
not a simple power law velocity dependence. Moreover, 
for the number of modes taken into account (n max = 40) 
the numerical results are not stable enough to resolve 
the weak dependence of m c on Vb with a high enough 
accuracy. (But it is good enough to perfectly resolve 
the ^-dependence.) The reason for the slow convergence 
of the numerics will become clear below. As we shall 
see, the corresponding energy density is dominated by 
masses much larger than m c . Consequently the weak 
dependence of m c on Vb is not very important in that 
respect and therefore does not need to be determined 
more precisely. However, combining all the data we can 
give as a fair approximation 



m c (k,v b ) 



with 



0.1. 



(5.3) 



Taking a = 0.13 for 0.1 < v b < 0.3 and a = 0.08 for 
0.03 < Vb < 0.1 fits the given data reasonably well. 

As we have seen, as long as the zero mode is the 
dominant source of KK-particle production, the final 
KK-graviton spectrum can be approximated by a power 
law decrease rn~ 3 . We can combine the presented 
numerical results to obtain a fitting formula valid in this 
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regime: 



a rout 

■ /v n»l,fc<l^ 



K) 



2.37 



k y s (L m ri 



for — < m r , 



(5.4) 



This fitting formula is shown in Figs. |2lj |22| and |23j 
and is in reasonabl e go od agreement w ith the numerical 
results. Since Eq. (5^) together with (5^3) is an impor- 
tant result, we have reintroduced dimensions, i.e. the 
AdS scale L which is set to one in the simulations, in 
both expressions. 

Let us now investigate the slope of the KK-graviton 
spectrum for masses m n — > oo since it determines the 
contribution of the heavy KK modes to the energy den- 
sity. In Figure ^5] we show KK-graviton spectra obtained 
for three- momentum k = 0.1, second brane position 
y s = 3 and maximal brane velocities Vb = 0.01,0.03 
and 0.1. Up to n max = 100 KK modes have been taken 
into account in the simulations. One immediately is 
confronted with the observation that the convergence 
of the KK-graviton spectra for large m n is very slow. 
This is since those modes, which are decoupled from the 
evolution of the four-dimensional graviton, are strongly 
affected by the kink in the brane motion. Recall that the 
production of light KK gravitons with masses m n <C m c 
is virtually driven entirely by the evolution of the 
massless mode. Those light modes are not so sensitive 
to the discontinuity in the velocity of the brane motion. 
To be more precise, their primary source of excitation is 
the evolution of the four-dimensional graviton but not 
the kink which, as we shall discuss now, is responsible 
for the production of heavy KK gravitons m n 3> m c . 
A discontinuity in the velocity will always lead to a di- 
vergent total particle number. Arbitrary high frequency 
modes are excited by the kink since the acceleration 
diverges there. Due to the excitation of KK gravitons 
of arbitrarily high masses, one cannot expect that the 
numerical simulations show a satisfactory convergence 
behavior which allows to determine the slope by fit- 
ting the data. However, it is nevertheless possible to 
give a quantitative expression for the behavior of the 
KK-graviton spectrum for large masses. The studies of 
the usual dynamical Casimir effect on a time-dependent 
interval are very useful for this purpose. 

For the usual dynamical Casimir effect it has been 
shown analytically that a discontinuity in the velocity 
will lead to a divergent particle number ]57], [5£|. In 
Appendix |] we discuss in detail the model of a massless 
real scalar field on a time-dependent interval [0,y(r,)] for 
the boundary motion y(t) = yo + vt with v = const, 
and present numerical results for final particle spectra 
(Fig. 34). For this motion it was shown in |58| that the 
particle spectrum behaves as oc v 2 /uj n where tu n = nir/yo 
is the frequency of a massless scalar particle. This di- 
vergent behavior is due to the discontinuities in the 
velocity when the motion is switched on and off, and are 
responsible for the slow convergence of the numerical 
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FIG. 25: KK-particle spectra for k = 0.1, y s = 3 and max- 
imal brane velocities vt = 0.01,0.03,0.1 up to KK masses 
m n ~ 100 compared with an l/m n decline. The dashed lines 



indicate the approximate expression (5.6) which describes the 
asymptotic behavior of the final KK-particle spectra reason- 
ably well, in particular for Vb < 0.1. 



results shown in Fig. [34| for this scenario. 
At the kink in the brane-motion the total change of the 
velocity is 2vb, similar to the case for the linear motion 
where the discontinuous change of the velocity is 2v. 
Consequently we may conclude that for large KK masses 
m n S> m c for which the evolution of the KK modes is no 
longer affected by their coupling to the four-dimensional 
graviton the KK-graviton spectrum behaves as 10 
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OC 



for m n 3> m c . 



(5.5) 



If we assume that the spectrum declines like l/m n and 
use that the numerical results for masses m n ~ 20 are 
virtually stable one finds Af^k • ^ v b'° 8 / m n which de- 
scribes the asymptotics of the numerical results well. 
As for the dynamical Casimir effect for a uniform motion 
discussed in Appendix |] [cf. Fig. [34j, the slow conver- 
gence of the numerical results towards the l/m n behav- 
ior is well visible for large masses m n ^> m c which do no 
longer couple to the four-dimensional graviton. This is a 
strong indication for the statement that the final gravi- 
ton spectrum for large masses behaves indeed like (j5~5|). 
It is therefore possible to give a single simple expression 
for the final KK-particle spectrum for large masses which 



10 Note that the discussion in Appendix |E] refers to Dirichlet bound- 
ary conditions. For Neumann boundary conditions considered 
here, the zero mode and its asymmetric coupling play certainly a 
particular role. However, as we have shown, for large masses only 
the KK-intermode couplings are important. Consequently, there 
is no reason to expect that the qualitative behavior of the spec- 
trum for large masses depends on the particular kind of boundary 
condition. 
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FIG. 26: Evolution of the zero-mode particle number A/"o,fc,. (t) 
and final KK-graviton spectra • f° r J/s = 3, maximal 
brane velocity Vb = 0.1 and three- momenta fc = 10 and 30. 



The dashed line in the upper plots indicate Eq. (5.17) (divided 
by two) demonstrating the value of the number of produced 
zero-mode gravitons without coupling to the KK modes. 



FIG. 27: 4D-graviton number Af ,k,.(t) for k = 3, 5, 10, 20 and 
30 with y s — 3 and maximal brane velocity Vb = 0.1. The 
small plot shows the final graviton spectrum A/Jfe*. together 
with a fit to the inverse law a l(k + b) [dashed line] and the 
analytical fitting formula Eq. (p.2S|) [solid line]. For k = 10 
and 30 the corresponding KK-graviton spectra are shown in 
Fig. H 



comprises all the features of the spectrum even quantita- 
tively reasonably well [cf. dashed lines in Fig. [25| 

Af™l. ~ 0.2^|— for m n » m c . (5.6) 

The l/?/ s -dependence is compelling. It follows imme- 
diately from the considerations on the energy and the 
scaling behavior discussed above [cf. Figs. [l6| and p"7| . 
For completeness we now write l/w° u ^ instead of the KK 
mass m n only, since what matters is the total energy 
of a mode. Throughout this section this has not been 
important since we considered only k <C 1 such that 
w™J becomes independent of k for large masses m n ^> k 
[cf.' Fig. [9|. 

E. Short wavelengths k ;§> 1 

For short wave lengths k ^> 1 (short compared to the 
AdS-curvature scale L set to one in the simulations) a 
completely new and very interesting effect appears. The 
behavior of the four-dimensional graviton mode changes 
drastically. We find that the zero mode now couples 
to the KK gravitons and no longer evolves virtually 
independently of the KK modes, in contrast to the 
behavior for long wavelengths. 

In Fig. ^6] we show the evolution of the zero-mode 
graviton number A/"o,fe,»(t) and final KK-graviton spectra 
■Af° u £. for y s = 3, maximal brane velocity v b = 0.1 
and three-momenta k = 10 and 30. One observes that 
the evolution of the four-dimensional graviton depends 
on the number of KK modes n max taken into account, 
i.e. the zero mode couples to the KK gravitons. For 



k = 10 the first 60 KK modes have to be included in the 
simulation in order to obtain a numerically stable result 
for the zero mode. In the case of k = 30 one already 
needs n max — 100 in order to achieve numerical stability 
for the zero mode. 

Figure ^7] displays the time-evolution of the number 
of produced zero-mode gravitons J\fo,k,»(t) for y s = 3 
and Vb = 0.1. For large k the production of massless 
gravitons takes place only at the bounce since these 
short wavelength modes are sub-horizon right after 
the bounce. Corresponding KK-particle spectra for 
k = 10, 30 are depicted in Figs. ^6] and ^8|. The insert 
in Fig. |27] shows the resulting final four-dimensional 
graviton spectrum A/J^*., which is very well fitted 
by an inverse power law A/" ™*. = 0.02/(fc - 1.8) 11 . 
Consequently, for k 3> 1 the zero-mode particle number 
A/'gfc*. declines like 1/fc only, in contrast to the 1/k 2 
behavior found for k -C 1. 

The dependence of A/J^*. on the maximal brane 
velocity Vb also changes. In Fig. we show Ao,fc,.(i) 
together with the corresponding KK-graviton spectra for 
y s = 3, k — 5 and 10 in each case for different V},. Using 
"■max = 60 KK modes in the simulations guarantees 
numerical stability for the zero mode. The velocity 
dependence of A/J ^ . is not given by a simple power law 
as it is the case for k <C 1. This is not very surprising 



1 The momenta k = 5, 10, 20, 30 and 40 have been used to obtain 
the fit. Fitting the spectrum for k = 20, 30 and 40 to a power 
law gives A/J^*, <x fc" 11 . 
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FIG. 28: Zero-mode particle number JVb,fc,»(i) and corre- 
sponding final KK-particle spectra Af°^,m for y s = 3, k = 5, 10 
and different maximal brane velocities vi,. n max = 60 guaran- 
tees numerically stable solutions for the zero mode. 



FIG. 29: Zero-mode particle number A/"o,fc.(t) and correspond- 
ing KK-graviton spectra for k = 10, vt, = 0.1 and 2nd brane 
positions y s — 3 and 10. 



since now the zero mode couples strongly to the KK 
modes [cf. Fig. [2q] . For k = 10, for example, one finds 

yy-out 
■ /v 0,fc,» 



oc vi A if v b < 0.1. 



As in the long wavelengths case, the zero-mode 
particle number does not depend on the position of the 
static brane y s even though the zero mode now couples 
to the KK modes. This is demonstrated in Fig. p9| where 
the evolution of the zero-mode particle number A<o,fc,»(i) 
and the corresponding KK-graviton spectra with A; = 10, 
Uf, = 0.1 for the two values y s = 3 and 10 are shown. 
One needs n max = 60 for y s — 30 in order to obtain a 
stable result for the zero mode which is not sufficient in 
the case y s — 10. Only for n max ~ 120 the zero-mode 
solution approaches the stable result which is identical 
to the result obtained for y s = 3. 

What is important is not the number of the KK modes 
the four-dimensional graviton couples to, but rather 
a particular mass m zm ~ k. The zero mode couples 
to all KK modes of masses below m zm no matter how 
many KK modes are lighter. Recall that the value of y s 
just determines how many KK modes belong to a given 
mass interval Am since, roughly, m n ~ mr/y s . The 
KK-spectra for k > 1 show the same scaling behavior as 
demonstrated for long wavelengths in Figs. |l6| and 17. 



The production of four-dimensional gravitons of 
short wavelengths takes place on the expense of the 
KK modes. In Fig. 30 we show the numerical results 
for the final KK-particle spectra with vj, — 0.1, y s = 3 
and k = 3,5, 10 and 30 obtained for different coupling 
combinations. These spectra should be compared with 
those shown in Fig. ^ for the long wavelengths case. For 
k > 10 the number of the produced lightest KK gravi- 
tons is smaller in the full coupling case compared to the 
situation where only the KK-intermode coupling is taken 
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FIG. 30: Final KK-particle spectra Af°% m for v b = 0.1, y a = 3 
and k = 3, 5, 10 and 30 and different couplings. Circles cor- 
respond to the full coupling case, squares indicate the results 
if Mij = Mu — 0, i.e. no KK-intermode couplings and dia- 
monds correspond to Mio = 0, i.e. no coupling of KK modes 
to the zero mode. 



into account. In case k = 30, for instance, the numbers 
of produced gravitons for the first four KK modes are 
smaller for the full coupling case. This indicates that 
the lightest KK modes couple strongly to the zero mode. 
Their evolution is damped and graviton production in 
those modes is suppressed. The production of zero-mode 
gravitons on the other hand is enhanced compared to 
the long wavelengths case. For short wavelengths, the 
evolution of the KK modes therefore contributes to 
the production of zero-mode gravitons. This may be 
interpreted as creation of zero-mode gravitons out of 
KK-mode vacuum fluctuations. 

As in the long wavelengths case, the KK-particle 
spectrum becomes independent of k if m n S> k and 
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FIG. 31: Final KK-particle spectra Af°%. for v b = 0.1, 
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simple analytical expression ( |6.35| ) agrees quite well with the 
numerical results. 



FIG. 32: KK-particle spectrum for y s = 3, Vb = 0.1 and 
k = 0.1 for the bouncing as well as smooth motions with 
is = 0.005, 0.015, and 0.05 to demonstrate the influence of the 
bounce. n max = 60 KK modes have been taken into account 
in the simulations and the result for the kink motion is shown 
as well. 



F. A smooth transition 



the evolution of the KK modes is dominated by the 
KK-intermode coupling. This is visible in Fig. [3^ for 
k = 3 and 5. Also the bend in the spectrum when the 
KK-intermode coupling starts to dominate is observable. 
For k = 10 and 30 this regime with m n ^s> k is not 
reached. 

As we have shown before, in the regime m n 3> k 
the KK-particle spectrum behaves as l/w° u ^ which will 
dominate the energy density of produced KK gravitons. 
If 1 <C m n < k, however, the zero mode couples to 
the KK modes and the KK-graviton spectrum does 
not decay like l/w° u ^. This is demonstrated in Fig. |3l] 
where the number of produced final state gravitons 
AT° U £ m is plotted as function of their frequency w° u ^ for 
parameters Vb — 0.1, y s — 3 and k = 5, 10, 20, 30 and 40. 
While for k = 5 the KK-intermode coupling dominates 
for large masses [cf. Fig. |o| leading to a bending over 
in the spectrum and eventually to an l/w° u ^-decay, the 
spectra for k = 20, 30 and 40 show a different behavior. 
All the modes are still coupled to the zero mode leading 
to a power-law decrease oc l/(w° u J) a with a ~ 2. The 
case k = 10 corresponds to an intermediate regime. 
Also show n is the simple analytical expression given in 
Eq. ( 6.3 5| ) which describes the spectra reasonably well 
for large k (dashed line). 

The KK-particle spectra in the region 1 <C m„ ^ k will 
also contribute to energy density since the cutoff scale is 
the same for the integration over k and the summation 
over the KK-tower (see Section VI D below). 



Let us finally investigate how the KK- graviton spec- 
trum changes when the kink- motion (2.18) is replaced by 
the smooth motion (5.1). In Fig. ^we show the numeri- 
cal results for the final KK-graviton spectrum for y s = 3, 
Vb = 0.1 and k — 0.1 for the smooth motion (5.1) with 
t s = 0.05, 0.015 and 0.005. n max = 60 modes have been 
taken into account in the simulation and the results are 
comp ared to the spectrum obtained with the kink-motion 
( |2.18| ). The parameter t s defines the scale L s ~ 2t s at 
which the kink is smoothed, i.e. L s corresponds to the 
width of the transition from contraction to expansion. 
The numerical results reveal that KK gravitons of masses 
smaller than m s ~ 1/L S are not affected, but the pro- 
duction of KK particles of masses larger than m s is 
exponentially suppressed. This is in particular evident 
for t s — 0.05 where the particle spectrum for masses 
m n > 10 has been fitted to a exponential decrease. Going 
to smaller values of t s , the suppression of KK-mode pro- 
duction sets in for larger masses. For the example with 
t s = 0.005 the KK-particle spectrum is identical to the 
one obtained with the kink-motion within the depicted 
mass range. In this case the exponential suppression of 
particle production sets in only for masses m n > 100. 
Note that the exponential decay of the spectrum for the 
smooth transition from contraction to expansions also 
shows that no additional spurious effects due to the dis- 
continuities in the velocity when switching the brane dy- 
namics on and off occur. Consequently, t- m and t ou t are 
appropriately chosen. 
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VI. ANALYTICAL CALCULATIONS AND 
ESTIMATES 



is just a delta- function potential with "height" propor- 
tional to —2y/Tib/L 



A. The zero mode: long wavelengths k <C 1/L 



(6.10) 



The numerical simulations show that the evolution of 
the zero mode at large wavelengths is not affected by the 
KK modes. To find an analytical approximation to the 
numerical result for the zero mode, we neglect all the 
couplings of the KK modes to the zero mode by setting 
Mij = Vi, j and keeping Mqo only. Then only the evo- 



lution equation for e Q 
and reduces to 



(a) 



5q6 is important; it decouples 



wher e v b i s given in Eq (2.20). Equation (6.1) with poten- 
tial ( J6. L0| ) can be considered as a Schrodinger equation 
with (5-function potential. Its solution is a classical text- 
book problem. 

Since the approximated potential V van ishes for allt < 
one has, with the initial condition (6.3), 



e(t) 



-ikt 



t < o 



(6.11) 



e+[k 2 + V(t)]e = , 



with "potential" 



(6.1) Assuming continuity of e through t — and integrating 
the differential equation over a small interval t € [0~ , + ] 
around t = gives 



V = M t 



,J oo — M 00 



(6.2) 



The corresponding vacuum initial conditions are 
[cf. Eqs. (3.21), (3.22); here we do not consider the unim- 
portant phase] 



lim e = 1 , lim e = — ik. 

t — > — oo t — > — oo 



(6.3) 



A brief calculation using the expression for M 00 (cf . Ap- 
pendix H) leads to 



V = 



vl 



Vb_ + Vb m 

Vb 



2y s 2 



9 9 



II- 



vl 



n 2 i 



-vl 



vl 



n 



(6.4) 
(6.5) 



If one assumes that the static brane is much further away 
from the Cauchy horizon than the physical brane, y s 3> 
y b , it is simply 



V 



-H 2 -H , 



(6.6) 



and one recovers Eq. (2.50) 



For the particular scale factor ( 2.17 ) one obtains 

sgn(t) 



H = - = 

a 



n 



t\ + t b 

2S(t) 1 



and 



t b i\t\+t b ) 2 



such that 



H + H 



2S(t) 



(6.7) 
(6.8) 

(6.9) 



The (5-function in the last equation models the bounce. 
Without the bounce, i.e. for an eternally radiation dom- 
inated dynamics, one has V = and the evolution equa- 
tion for e would be trivial. With the bounce, the potential 



= 



fc 2 - ^ff(t) . < 



e(0 + )-e(0_)-^e(0) 



(6.12) 
(6.13) 



The jump of the derivative e at t = leads to parti- 
cle creation. Using e(0+) = e(0) = e(0_) and e(0 + ) = 
e(0_) + 2 -^e(0) as initial conditions for the solution for 
t > 0, one obtains 



e{t) = Ae~ lkt + Be 



i k t 



t > o 



with 



A= 1+i 



Vb 

kL 



B 



Vb 

kL 



(6.14) 



(6.15) 



The Bogoliubov coefficient Bqq after the bounce is then 
given by 



B Q0 {t > 0) 



-ikt 



(6.16) 



where we have used that Moo = ~ W if y s ^> yb- At 
this point the importance of the coupling matrix Moo 
becomes obvious. Even though the solution e to the dif- 
ferential equation (3.1) is a plane wave right after the 
bounce, |2?oo(£)| 2 is n °t a constant due to the motion of 
the brane itself. Only once the mode is inside the hori- 
zon, i.e. Ti/k -C 1, |$oo(i)| 2 i s constant and the number 
of generated final state gravitons (for both polarizations) 
is given by 



\f° l 



2\B 00 (kf> 1)| 2 

2v b 
(kL) 2 



k 2 



1 

2 

(6.17) 



where we have used that the Wronskia n of e, e* is 2ik. 
As illustrated in Fig. ||the expression ( |6.17 ) is indeed in 
excellent agreement with the (full) numerical results, not 
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only in its fc-dependence but also the amplitude agrees 
without any fudge factor. The evolution of the four- 
dimensional graviton mode and the associated genera- 
tion of massless gravitons with momentum k < 1/L can 
therefore be understood analytically. 
Note that the approximation employed here is only valid 
if y 2 — yb{0) 2 ^> Vb(0) 2 - In the opposite limit, if Ay = 
Vs — Ub{0) <C y&(0) one can also derive an analytical ap- 
proximation along the sam e lines. For k < 1/Ay one 
obtains instead of Eq. (3.17) 



For the ~ sign we assume t ^> L and t 3> 4 so that 
one may neglect terms of order t/L in comparison to 
^/v~b(t/L) 2 . We have also a ppr oximated a = (t + tb)/L ~ 
t/L. Inserting this in Eq. ( |4.§| ) yields 



V (k) 



«4 

2^3 



Vb k 



fcf < 1 



(6.22) 



Both expressions ( |6.2C ) and (6.22) are in very good 
agreement with the corresponding numerical results, see 
Figs, g, and [ll]. 



A/" 01 
■ /v 0,, 



2(kAy) 2 ' 
if Ay = y s - y b (0) < y b (0) 



(6.18) 
kAy < 1 . 



In order to calculate the energy density, we have to take 
into account that the approximation of an exactly ra- 
diation dominated Universe with an instant transition 
breaks down on small scales. We assume this break 
down to occ ur at t he str ing scale L s , much smaller than 
L [cf. Eqs. ([04l), (|2.15|) l. L s is the true width of the 
transition from collapse to expansion, which we have set 
to zero in the treatment. Modes with mode numbers 
k ^> (2tt)/L s will not 'feel' the potential and are not 
generated. We therefore choo se fc m ax = (2n)/L s as the 
cutoff scale. Then, with Eq (4.21), one obtains for the 
energy density 



1 



Po 



2 ir 2 a 4 



2k/L s 



dkk 3 J\f Q k 



(6.19) 



For small wave numbers, k < 1/L, we can use the above 
analytical result for the zero- mode particle number. 
However, as the numerical simulations have revealed, as 
soon as k > 1/L, the coupling of the four-dimensional 
graviton to the KK modes becomes important and for 
large wave numbers decays only like 1/k. Hence 



the integral (6.19) is entirely dominated by the upper 
cutoff. The contributions from long wavelengths to the 
energy density are negligible. 

For the power spectrum, on the other hand, we 
are interested in cosmologically large scales, 1/k ~ 
several Mpc or more, but not in short wavelengths 
kL S> 1 dominating the energy density. Inserting the 
expression for the number of produced long wavelength 
gravitons (6.17) into (4.11), the gravity wave power 
spectrum at late times becomes 



V {k) 



2 Vb K4 

W (aLf 



for kt > 1. 



(6.20) 



This is the asymptotic power spectrum, when e starts 
oscillating, hence inside the Hubble horizon, k > 1. 
On super Hubble scales, kt <§; 1 when the asymptotic 
out-state of t he zero mode is not yet reached, one may 
use Eq. ( gig ) with 



Av b a 2 



(6.21) 



B. The zero mode: short wavelengths 1/L 

As we have demonstrated with the numerical analysis, 
as soon as k > 1/L, the coupling of the zero mode to 
the KK modes becomes important, and for large wave 
numbers Af§^, cx 1/k. We obtain a good asymptotic 
behavior for the four-dimensional graviton spectrum if 
we set 
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5(kL) 



(6.23) 



This function and Eq. ( 6.17 ) (divided by two for one po- 
larization) meet at kL = 5. Even though the approxi- 
mation is not good in the intermediate regime it is very 
reasonable for large k [cf. Fig. E7|. 



Inserting this approximation into Eq (6.19) for the energy 
density, one finds that the integral is dominated entirely 
by the upper cutoff, i.e. by the blue, high energy modes: 



16 7T Vb 1 IT Vb 



30 a 4 LL\ 2 a 4 LL\ ' 



(6.24) 



The power spectrum associated with the short wave- 
lengths k ^> 1/L is not of interest since the gravity wave 
spectrum is measured on cosmologically large scales only, 
k < 1/L. 



C. Light Kaluza-Klein modes and long wavelengths 

fc<l/L 

The numerics indicates that light (m„ < 1) long wave- 
length KK modes become excited mainly due to their 
coupling to the zero mode. Let us take only this cou- 
pling into account and neglect also the time-dependence 
of the frequency, setting to n ,k{t) = w™t = tu™ k since it 
plays an inferior role as shown by the numerics. 
The Bogoliubov coefficients are then determined by the 
equations 



2w ou l 

n.k 



71, K 



S n (t;k) (6.25) 



S n {t;k) (6.26) 



with the "source" 



S n (t;k) = (£ -i]o)M n0 



(6.27) 
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We have defined £„ 



AO) „ , = „(o) t = *(o) „ nH 

?n,fe' Vn.k — V n ,k> SO — ? ,fc' ailCL 



r/o = ^nfe" This source is known, since the evolution of 
the four-dimensional graviton is know. From the result 
for e above and the definition of £o and ?7o in terms of e 
and e one obtains 



c 21 



1 



t\+t b 



-itk 



, t < (6.28) 
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1 — ifcth 



kt b k 2 t b (t + t b ) 
i 1 



-ilk 



kt b kH b (t + t b ) 



itk 



t > (.6.29) 



Furthermore, if y s ^> y bl one has [cf. Eq. (B3)] 



M, 



Y 1 (m n y s ) 2 

Vb V Yi(m n y b ) 2 - Yi(m n y s y 



(6.30) 



Assuming y s m n 3> 1 and y b m n <C 1 one can expand the 
Bessel functions and arrives at 




' irm n L 2 Lsgn(t) 

y s (\t\ + h) 2 



To determine the number of created final s tate gravi- 
tons we only need to calculate rj n f. [cf. Eq. ( 3.32 ) with 

- °°) = 1 and A nA\t\ -'«>) = 0], 



■ACfe,. = |^0n,fc(iout)| 2 



iw: 



4 & 



Vn,k\ 



(6.31) 



The vacuum initial conditions require lim^-oo r\ n ^ = 
so that rj n /. is given by the particular solution 



J n,k 



5 i (^ , ;/fc)e- ^t '^c^^ , 



and therefore 



4<£ 



S n {t;k)e- itu]a ^dt 



(6.32) 



(6.33) 



where the integration range has been extended from — oo 
to +oo since the source is very localized around the 
bounce. This integral can be solved exactly. A some- 
what lengthy but straight forward calculation gives 



ACL = 



^5 t4 



2iRc (e'^+^^^w^ + k)t b ) 



+(tt i )- 1 e'«-H^ 1 ( ! ( w - _ k )t b ) 



-e t( -<" +k ^E 1 (i(u J n ut k + k)t b ) 



(6.34) 



Here E\ is the exponential integral, E\(z) = 
joo f.-ig-tfj^ _ This function is holomorphic in the com- 
plex plane with a cut along the negative real axis, and 
the above expression is therefore well defined. Note that 



this expression does not give rise to a simple dependence 
of Af° u £ on the velocity v b = (L/t b ) 2 . In the preceding 
secti on w e have seen that, within its range of validity, 



Eq. ( 3.34 ) is in excellent agreement with the numerical 
results (cf., for instance, Figs. [l2| and p^| ). 
As already mentioned before, this excellent agreement 
between the numerics and the analytical approximation 
demonstrates that the numerical results are not contam- 
inated by any spurious effects. 



D. Kaluza-Klein modes: asymptotic behavior and 
energy density 

The numerical simulations show that the asymptotic 
KK-graviton spectra (i.e. for masses m n ^> 1) decay 



with a 



like l/tJ° w l if m n > k and like 
if m n < k. The corresponding energy density on the 
brane is given by the summation of Eq. (4.23) over all 
KK modes up to the cutoff. Since the mass m n is simply 
the momentum into the extra dimension, it is plausible 
to choose the same cutoff scale for both, the fc-integral 
and the summation over the KK modes, namely 2tt/L s . 
The main contribution to the four-dimensional particle 
density and energy density comes from m n ~ 2n/L s and 
k ~ 2ir/L s , i.e. the blue end of the spectrum. 
The large-frequency behavior of the final KK-spectrum 
can be approximated by 



« rout 

• /v n,/c,. 



0.2v 2 



1 



J n,k 



if l/L<k< m r . 



-l)/2 



k a ~ 



if 



< k < 2tt/L s 



(6.35) 

with a ~ 2 which is particularly good for large k. Both 

fression match at m n = k and are indicated in Figures 
and " 
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as dashed lines. Given the complicated 
coupling structure of the problem and the multitude of 
features visible in the particle spectra, these compact 
expressions describe the numerical results reasonable 
well for all parameters. The deviation from the numer- 
ical results is at most a factor of two. This accuracy 
is sufficient in order to obtain a useful expression for 
the energy density from which bounds on the involved 
energy scales can be derived. 

The energy density on the brane ass ociated with 
the KK gravitons is given by [cf. Eq. ( 4.23| )1 



Pkk 



ira 6 y s 



E 



dkk 2 A/? ut 



" n,k,t 



u n,k 



(6.36) 



Splitting the momentum integration into two integrations 
from to m n and m n to the cutoff 2n/L s , and replacing 
the sum over the KK masses by an integral one obtains 



.32 



12 



Pkk - C(a) 



a 6 y s L 5 S 



(6.37) 



The power a in Eq. (5.35) enters the final result for the 
energy density only through the pre-factor C (a) which is 
of order unity. 



VII. DISCUSSION 

The numerical simulations have revealed many inter- 
esting effects related to the interplay between the evolu- 
tion of the four-dimensional graviton and the KK modes. 
All features observed in the numerical results have been 
interpreted entirely on physical grounds and many of 
them are supported by analytical calculations and ar- 
guments. Having summarized the results for the power 
spectrum and energy densities in the preceding section, 
we are now in the position to discuss the significance of 
these findings for brane cosmology. 



A. The zero mode 



For the zero-mode power-spectrum we have found that 
'- 2 if kt < 1 



(7.1) 



Therefore, the gravity wave spectrum on large, super 
Hubble scales is blue with spectral tilt 



TIT 



(7.2) 



a common feature of ekpyrotic and pre-big-bang models. 
The amplitude of perturbations on scales at which 
fluctuations of the Cosmic Microwave Background 
(CMB) are observed is of the order of (H /mpi) 2 , i.e. 
very suppressed on scales relevant for the anisotropics 
of the CMB. The fluctuations induced by these Casimir 
gravitons are much too small to leave any observable 
imprint on the CMB. 

For the zero-mode energy de nsity at late times, 
kt » 1, we have obtained [cf Eq. ( 5.24 )1 



1 7T Wfc 

ph °- 2^zzr 



(7.3) 



In this section we denote the energy density of the 
zero mode by p^o in order not to confuse it with the 



12 Note that even the transition from the sum mation over the KK- 
tower to an integration according to ( k-.3,3TI "eats up" the l/y s 
term in (6.3C ), the final energy density (|6.37|) depends on y a since 



present density of the Universe. Recall t hat L s is the 



scale at which our kinky approximation ( 2.17 ) of the 
scale factor breaks down, i.e. the width of the bounce. 
If this width is taken to zero, the energy density of 
gravitons is very blue and diverges. This is not so 
surprising, since the kink in a(t) leads to the generation 
of gravitons of arbitrary high energies. However, as the 
numerical simulations have shown, when we smooth the 
kink at some scale L s , the production of modes with 
energies larger than ~ 1/L S is exponentially suppressed 
[cf. Fig. p2[ . This justifies the introduction of L s as a 
cutoff scale. 

In the following we shall determine the density pa- 
rameter of the generated gravitons today and compare 
it to the Nucleosynthesis boun d. For this we need the 
quantities at, given in Eq ( 2.2C ) and 



H h = 



Vb 

L 



Here Of, is the minimal scale factor and is the maximal 
Hubble parameter, i.e. the Hubble parameter right after 
the bounce. (Recall that in the low energy approximation 
t = r/.) During the radiation era, curvature and/or a cos- 
mological constant can be neglected so that the density 
is 



Prad 



K-4 



K4L 2 



h'4 



-Hi 



(7.4) 



In order to determine the density parameter of the gen- 
erated gravitons today, i.e., at t = to, we use 



n 



hO 



Pho(to) Pho{to) Prad(to) Pho(to) 



Pcrit(^o) Prad(ifj) Pcrit(^o) Prad(^o) 



rad • 



(7.5) 

The second factor f2 ra d is the present radiation density 
parameter. For the factor pho/Pr&d, which is time inde- 
pendent since both pho and p ra d scale like 1/a 4 , we insert 
the above results and obtain 



y ^rad 



n PhO n ITT ( Lpi 

"/i0 — "rad — n~o Vb 1 — 

Prad ^ O \ -Ls 



The nucleosynthesis bound [|14) requests that 

&h0 ^3 0.1 f2 ra d , 

which translates into the relation 

^(L Pl /L s f(L/L s ) <0.1 



(7.6) 
(7.7) 

(7.8) 
(7.9) 



which, at first sight, relates the different scales involved. 
But since we have chosen the cutoff scale L s to be 
the higher -dim ensional fundamental scale (string scale), 
Equation (7.9) reduces to 



it explicitly enters the particle number. 



v b < 0.2 



(7.10) 
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by virtue of Equation (2.15). All one has to require to 
be consistent with the nucleosynthesis bound is a small 
brane velocity which justifies the low energy approach. 
In all, we conclude that the model is not severely con- 
strained by the zero mode. This result itself is remark- 
able. If there would be no coupling of the zero mode 
to the KK modes for small wavelengths the number of 
produced high energy zero-mode gravitons would behave 
as oc fc~ 2 as it is the case for long wavelengths. The 
production of high energy zero-mode gravitons from KK 
gravitons enhances the total energy density by a factor of 
about L/L s . Without this enhancement, the nucleosyn- 
thesis bound would not lead to any meaningful constraint 
and would not even require Vb < 1. 



B. The KK modes 

As derived above, the energy density of KK gravitons 
on the brane is dominated by the hig h ene rgy gravitons 
and can be approximated by [cf. Eq. (6.37)] 



Pkk 



n 5 v 2 L 2 
a 6 y s L 5 S 



(7.11) 



Let us evaluate the constraint induced from the require- 
ment that the KK-energy density on the brane be smaller 
than the radiation density pkk(^) < Prad(i) at all times. 
If this is not satisfied, back-reaction cannot be neglected 
and our results are no longer valid. Clearly, at early times 
this condition is more stringent than at late times since 
Pkk decays faster then p ra d- Inserting the value of the 
scale factor directly after the bounce where the produc- 



tion of KK gravitons takes place, a b i 
using again the RS fine tuning condition (2.15) 



Vb, 



one finds, 



Pkk 

Prad 



lOOvi I — I I — 

Vs 



(7.12) 



If we use the largest value for the brane velocity Vb ad- 
mitted by the nucleosynthesis bound Vb — 0.2 and re- 
quire that Pkk/ Prad be (much) smaller than one for back- 
reaction effects to be negligible, we obtain the very strin- 
gent condition 



L L s 
— « — 

Vs \L 



(7.13) 



Let us first discuss the largest allowe d va lue for L ~ 
0.1mm. The RS-fine tuning condition (2.15) then deter- 



mines L s = (LLp;) 1 / 3 ~ 10~ 22 mm ~ 1/(10 6 TeV). In 



this case the brane tension is T = 6K4/K 2 = QL 2 pl /L & s = 
&/(LL 3 s ) ~ (10 TeV) 4 . Furthermore, we have (L/L s ) 2 ~ 
10 42 so that y 8 > L{L/L S ) 2 ~ 10 41 mm ~ 3 x 10 15 Mpc, 
which is about 12 orders of magnitude larger than the 
present Hubble scale. Also, since y b (t) <C L in the low 
energ y regime, and y s 3> L according to the inequality 
( [7.13 ), the physical brane and the static brane are very 



far apart at all times. Note that the distance between 
the physical and the static brane is 



-dy = L\og(y s /y b ) > L > L s 



y 



This situation is probably not very realistic. Some high 
energy, stringy effects are needed to provoke the bounce 
and one expects these to be relevant only when the 
branes are sufficiently close, i.e. a t a d istance of order 
L s . But in this case the constraint (7.13) will be violated 
which implies that back-reaction will be relevant. 

On the other hand, if one wants that ys ^ L and 

back-reaction to be unimportant, then Eq. ( 7.12D implies 
that the bounce velocity has to be exceedingly small, 
v b < 10~ 15 . 

A way out of this conclusion is to assume that the 
brane distance at the bounce, Ay = y s — yb(0), becomes 
of the order of the cutoff L s or smaller. Then the pro- 
duction of KK grav itons is suppressed. However, then 
the approximation (6.18) has to be used to determine 
the energy density of zero-mode gravitons which then 
becomes 

p m ^ V l{L s Ay)-\ 

Setting Ay ~ L s , the nucleosynthesis bound, 
Pho ~ 0.1p ra d, then yields the much more stringent limit 
on the brane velocity, 



L 



(7.14) 



One might hope to find a way out of these conclusions 
by allowing the bounce to happen in the high energy 
regime. But then Vb — 1 and the nucleosynthesis bound 
is violated since too many zero-mode gravitons are pro- 
duced. Even if one disregards this limit for a moment, 
saying that the calculation presented here only applies in 
the low energy regime, v b -C 1, the modification coming 
from the high energy regime are not expected to allevi- 
ate the bounds. In the high energy regime one may of 
course have y b (t) 3> L and therefore the physical brane 
can approach the static bra ne arb itrarily closely without 
the latter having to violate ( 7.13 ) . Those results suggest 
that even in the scenario of a bounce at low energies, the 
back reaction from KK gravitons has to be taken into 
account. But this docs not need to exclude the model. 



VIII. CONCLUSIONS 

We have studied the evolution of tensor perturbations 
in braneworld cosmology using the techniques developed 
for the standard dynamical Casimir effect. A model 
consisting of a moving and a fixed 3-brane embedded in 
a five-dimensional static AdS bulk has been considered. 
Applying the dynamical Casimir effect formulation 
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to the study of tensor perturbations in braneworld 
cosmology represents an interesting alternative to 
other approaches existing in the literature so far and 
provides a new perspective on the problem. The explicit 
use of coupling matrices allows us to obtain detailed 
information about the effects of the intcrmodc couplings 
generated by the time-dependent boundary conditions, 
i.e. the brane motion. 

Based on the expansion of the tensor perturbations 
in instantaneous eigenfunctions, we have introduced a 
consistent quantum mechanical formulation of graviton 
production by a moving brane. Observable quantities 
like the power spectrum and energy density can be 
directly deduced from quantum mechanical expectation 
values, in particular the number of gravitons created 
from vacuum fluctuations. The most surprising and 
at the same time most interesting fact which this 
approach has revealed is that the energy density of the 
massive gravitons decays like 1/a 6 with the expansion 
of the Universe. This is a direct consequence of the 
localization of gravity: five-dimensional aspects of it, 
like the KK gravitons, become less and less 'visible' 
on the brane with the expansion of the Universe. The 
l/a 6 -scaling behavior remains valid also when the fixed 
brane is sent off to infinity and one ends up with 
a single braneworld in AdS, like in the original RS II 
scenario. Consequently, KK gravitons on a brane moving 
through an AdS bulk cannot play the role of dark matter. 

As an explicit example, we have studied graviton 
production in a generic, ekpyrotic-inspired model of 
two branes bouncing at low energies, assuming that 
the energy density on the moving brane is dominated 
by a radiation component. The numerical results have 
revealed a multitude of interesting effects. 
For long wavelengths kL <C 1 the zero mode evolves 
virtually independently of the KK modes, zero-mode 
gravitons are generated by the self coupling of the zero 
mode to the moving brane. For the number of produced 
massless gravitons we have found the simple analytical 
expression 2vb/(kL). These long wavelength modes 
are the once which are of interest for the gravitational 
wave power spectrum. As one expects for an ekpyrotic 
scenario, the power spectrum is blue on super-horizon 
scales with spectral tilt jit — 2. Hence, the spectrum 
of these Casimir gravitons has much too little power 
on large scales to affect the fluctuations of the cosmic 
microwave background. 

The situation changes completely for short wavelengths 
kL 3> 1. In this wavelength range, the evolution of the 
zero mode couples strongly to the KK modes. Produc- 
tion of zero-mode gravitons takes place on the expense of 
KK-graviton production. The numerical simulation have 
revealed that the number of produced short-wavelength 
massless gravitons is given by 2vb/(5kL). It decays 
only like 1/fc instead of the l//c 2 -behavior found for 
long wavelengths. These short wavelength gravitons 
dominate the energy density. Comparing the energy 



density with the nucleosynthesis bound and taking the 
cutoff scale to be the string scale L s , we have shown that 
the model is not constrained by the zero mode. As long 
as Vb < 0.2, i.e. a low energy bounce, the nucleosynthesis 
bound is not violated. 

More stringent bounds on the model come from 
the KK modes. Their energy density is dominated by 
the high energy modes which are produced due to the 
kink which models the transition from contraction to 
expansion. Imposing the reasonable requirement that 
the energy density of the KK modes on the brane be 
(much) smaller than the radiation density at all times in 
order for back reaction effects to be negligible, has led to 
two cases. On the one hand, allowing the largest values 
for the AdS curvature scale L ~ 0.1mm and the bounce 
velocity Vb — 0.2, back reaction can only be neglected 
if the fixed brane is very far away from the physical 
brane y s ~ 10 41 mm. As we have argued, this is not very 
realistic since some high energy, stringy effects provoking 
the bounce are expected to be relevant only when the 
branes are sufficiently close, i.e. y s ~ L s . On the other 
hand, by only requiring that y s ~ L 3> L s , the bounce 
velocity has already to be exceedingly small, Vb 10~ 15 , 
for back reaction to be unimportant. Therefore, one of 
the main conclusions to take away from this work is that 
back reaction of massive gravitons has to be taken into 
account for a realistic bounce. 

Many of the results presented here are based on 
numerical calculations. However, since the used ap- 
proach provides the possibility to artificially switch 
on and off the mode couplings, we were able identify 
the primary sources driving the time evolution of the 
perturbations in different wavelength and KK mass 
ranges. This has allowed us to understand many of the 
features observed in the numerical results on analytical 
grounds. 

On the other hand, it is fair to say that most of the 
presented results rely on the low energy approach, i.e. on 
the approximation of the junction condition (generalized 
Neumann boundary condition) by a Neumann boundary 
condition. Even though we have given arguments for 
the goodness of this approximation, it has eventually 
to be confirmed by calculations which take the exact 
boundary condition into account. This is the subject of 
future work. 



Acknowledgment 

We thank Cyril Cartier who participated in the early 
stages of this work and Kazuya Koyama and David Lan- 
glois for discussions. We are grateful for the use of the 
'Myrinet '-cluster of Geneva University on which most of 
the quite intensive numerical computations have been 
performed. This work is supported by the Swiss National 
Science Foundation. 



35 



APPENDIX A: VARIATION OF THE ACTION 



with 



Let us consider the variation of the action ( |2.27 ) with 
respect to h m . It is sufficient to study the action for a 
fixed wave number k and polarization • 



dt 



" dy 

,(*) y 3 



[\d t h.\ 2 -\d y h.\ 2 -k 2 \h.\ 2 ] 

(Al) 



and we omit the normalization factor L 3 //eg as well as 
the factor two related to Z 2 symmetry. The variation of 
(Al) reads 



SS h .(k) = U '*[" % 

1 JT Jy b {t) y 



(dth.)(d t 6K) 



-(dyh.)(d y Shl) - k 2 h.5hl 



(A2) 



h.c. 



Here, T denotes a time interval within the variation is 
performed and it is assumed in the following that the 
variation vanishes at the boundaries of the time interval 
T. Performing partial integrations and demanding that 
the variation of the action vanishes leads to 



= 



dt 



>> dy 

.(*) y 3 



{-d 2 h. 



h. 



y \y 3 



k 2 h. 



(A3) 



dt{± [{vdt + d y )h.]SK\ Vb{t) - ^(d y h.)Sh:\ ys } 



with v — dy b (t) /dt. T he first term in curly brackets is the 
wave operator ( ]2.24j) . In order for h, to satis fy the free 
wave equation (perturbation equation) ( 2.24 ) the term 
in curly brackets in the second integral has to vanish. 
Allowing for an evolution of h, on the branes, i.e. in 
general <S/i«|brane 7^ 0, enforces the boundary conditions 



(vd t + d y ) h.\ Vb{t) = and d y h.\ Vs = , (A4) 

hence, the junction condition ( ^.2(j) Consequently, any 
other boundary conditions than (A4) a rc no t compati- 
ble with the free perturbation equation ( 2.24 ) under the 
influence of a moving brane (provided 8h t ^ at the 
branes). 



APPENDIX B: COUPLING MATRICES 



The use of several identities of Bessel functions leads 



to 



M m = y b s „ , 

vi - vi 

M 0j = , 

4A^i i)b 4 y s 
M i0 = 4> = y b N, 



mm y b 
M u = fij , 
M i:j = M-j + 



(Bl) 
(B2) 
(B3) 

(B4) 
(B5) 



= [y b + m t )y b 



2m 2 NiNj 

m - — mf 

j 1 



(B6) 



x [y s C 2 (m j y s ) Jiim.ys) - y b C 2 {m :j y b ) J-iirrHVb)] 
where 



Ji(m 4 y) = [J2(rriiy b )Yi(m t y) - Y 2 {m l y b )Ji{m l y)} 

(B7) 



and 



M^NiNjimm d»irCi(mi»)Ca(mjy). (B8) 
Jyb 

This integral has to be solved numerically. Note that, 
because of the boundary conditions, one has the identity 

dyy 2 Ci{m l y)C 2 {m j y) = - I dyy 2 Ci{m l y)C {m j y). 
'yb Jyb 

Furthermore, one can simplify 



(B9) 



rr 1 s 2 2 Y x {m iVs ) 

Ji{m l y b ) = , Ji(m l y s ) = — 

nm.yt, nm.yh Yi(miy b ) 

(BIO) 

where the limiting value has to be taken for the last term 
whenever Yi(miy b ) = Y 1 (m l y s ) = 0. 



APPENDIX C: ON POWER SPECTRUM AND 
ENERGY DENSITY CALCULATION 

1. Instantaneous vacuum 



In Section III the in - out state approach to particle 
creation has been presented. The definitions of the in 
- and out- vacuum states Eq. (3.9) are unique and the 
particle concept is well defined and meaningful. 



If we interpret t ou t as a continuous time varia ble t , 
we can write the Bogoliubov transformation Eq. (3.24) 
as 





(CI) 

where at any time we have introduced a set of operators 
{»Q;,k»(i), o„ k ,(£)}■ Vacuum states defined at any time 
can be associated with these operators via 



aa, k ,.(*)|0,t) =0 Va,k 



(C2) 



Similar to Eq. ( 3.11| ) a "particle number" can be intro- 
duced through 



A/;,fc(<) = ^(0,m|a^.(*)a a , k ,.(*)|0,in) 

= 2]T|ZW(i)| 2 . 





(C3) 



3G 



We shall denote \0,t) as the instantaneous vacuum state 
and the quantity J\f a ,k(~t) as instantaneous particle num- 
ber 13 . However, even if we call it "particle number" and 
plot it in section ^ for illustrative reasons, we consider 
only the particle definitions for the initial and final state 
(asymptotic regions) outlined in section [II as physically 
meaningful. 



2. Power spectrum 



In order to calculate the power spectrum Eq. (4.7) we 
need to evaluate the expectation value 



(h m (t,y b ,k)hl(t,y b ,k')) ia = 

L 3 



(C4) 



where 
(0,in|.. 



introduced the shortcut 



we have 

.|0, in). Using the expansion ( 3.15 ) of q a ',w,»{t) 
in initial state operators and complex functions elfc(i) 
one finds 



(fe,k,.(t)C'k'.(*))i 



E 





2w 



<5^(k-k') 



0,k 



(C5) 

From the initial conditions (3.21) it follows that the sum 
in (C4) diverges at t = t in . This divergence is related to 
the usual normal ordering problem and can be removed 
by a subtraction scheme. However, in order to obtain a 
well defined power spectrum at all times, it is not suffi- 
cient just to subtract the term (l/2)(5 aa '/cu i ^ k )5^(k — 

k') which corresponds to (g a ,k,.(*m)<Za,/ k ' .(*in)) 

i n in the 

above expression. In order to identify all terms contained 
in the power spectrum we use the instantaneous particle 
concep t wh ich all ows u s to treat the Bogoliubov coeffi- 
cients (3.25) and ( [3.26 ) as continuous functions of time. 



First we express the complex functions e^ k in (C5) in 
terms of A^ a ,k{t) and B ia ^{t). This is of course equiv- 
alent to calculating the expectation value ( |C5| ) using 
[cf. Eq.@] 



<la,k,»(t) — 



1 



+ <-k,.(*)e«,*(*)l ( C6 ) 



and the Bogoliubov transformation Eq. JC1| ) . The result 
consists of terms involving the Bogoliubov coefficients 
and the factor (l/2)(S aa > /ui atk (t))8^ - k'), leading 
potentially to a divergence at all times. This term cor- 
responds to (0, i|<7a,k,.(i)<v k' . (^)|0, t), and is related to 



the normal ordering problem (zero-point energy) with re- 
spect to the instantaneous vacuum state |0, t). It can be 
removed by the subtraction scheme 

in,phys 

(C7) 

= ($«,k,.(f)& )k> (t)>in- (0,t|fe,k,.(*)6,k',.(*)|0,*) 

where we use the subscript "phys" to denote the physi- 
cally meaningful expectation value. 



Inser ting this expectation value into (C4), and using 
Eq. (O, we find 



(/i.(i,yi>,k)M£,yf,,k'))i n 
ar L L — ' 



(C8) 



with TZ a>k (t) defined in Eq. (p~9[). The function Q^ k 
appearing in Eq. ([4.9[) is explicitely given by 



C% h = 2 5R {&a,*Ap°,kBfia,k + ®«,k £ x 

'w fl ,t y a > (a) 



Wa'.fc y a {a) 



[ a',fc^/3a^/3Q' + Q a' ,k-A/3 a Bp a >] } 

(C9) 



and k appearing in Eq. (4.10) reads 



(CIO) 



3. Energy density 

In order to calculate the energy density we need to 

evaluat e the expectation value (hij(t, x, y b )h lJ (t, x, y b ))m- 
Using ( |2.22 ) and the relation e\A— k) = (ej-(k))* we ob- 
tain 

(/i y -(t,x,?/ fc )ft, y (i,x,y b )) in = ^2 ( - d k -<~ d ^ tt. > 



(2tt) 3 / 2 (2tt) 3 / 2 



(Cll) 



(ft.(t, Wj k)^(t, W> k / )) ine < ( k - k ')-X-(k) (e'' y '(k'))* 



By means of the expansion (3.17) the expectation value 
(h. (t, y b , k)h[, (t, y b , k')) in becomes 



(h.(t,y b ,k)hl{t,y b> k')) 
~ L 3 



(C12) 



13 It could be interpreted as the number of particles which would 
have been created if the motion of the boundary (the brane) 
stops at time t. 



From the definition of p a M,*{t) m Eq. (3.18) it is clear 
that this expectation value will in general contain terms 
proportional to the coupling matrix and its square when 
expressed in terms of e„ L However, we are interested in 
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the expectation value at late times only when the brane 
moves very slowly such that the mode couplings go to 
zero and a physical meaningful particle definition can be 
given. In this case we can set 



(p«,k,.(t)pLk',.'(*))i n = (fe,k,.(i)<7Lk',.'(*)). -(C13) 



Calculating this expectation value by using Eq. ( 3.15| ) 
leads to an expression which, as for the power spec- 
trum calculation before, has a divergent part related to 
the zero-point energy of the instantaneous vacuum state 
(normal ordering problem). We r emo ve this part by a 
subtraction scheme similar to Eq (|C7|). The final result 
reads 



in,phys 



e m (t)i w (t) 

1 = ^a.kV-jOac 

, ,in , ,in 

v/^.fe^.fc' 



(C14) 
5..,6^ (k-k<). 



Inserting this result into Eq. ( |CI2 ), splitting the summa- 
tions in sums over a = a' and a ^ a' and neglecting the 
oscillating a a' contributions (averaging over several 
oscillations), leads to 



(h.(t,y b ,k)kl(t,y b) k% 



(C15) 



1 Ks 



where the function K, a ^{t) is given by 



u /3,k 



(C16) 

and we have made use of Eq. (4.2). The relation be- 



tween J2p l^akW^/^dk anc ^ * ne num ber of created par- 
ticles can easily be established. Using this expression in 
Eq. ( |C11| ) leads eventually to 



(hij(t,yL,y b )h i:, (t,x,y b ))^ 



(C17) 



_1_ «5 \ ^ 
t2 T. 2-^ 



d 3 k 



■ic a , k {t)yl{a) 



a- L ^ J (2tt) 3 
where we have used that the polarization tensors satisfy 
]TeV(k) (e'«(k))*=2. (C18) 



The final expression for the energy density Eq. ( 4.1 q ) is 
then obtained by exploiting that n^/L = K4. 



APPENDIX D: NUMERICS 

In order to calculate the number of produced gravitons 
the system of coupled differential equations (3.34) and 



( |3.35 ) is solved numerically. The complex functions £ 



(/3) 

^cf! are decomposed into their real and imaginary parts: 

M _ U (P) . iv (P) JP) _ X V» , iv (P) (m) 

S Q ,fc — a a,k ^ lU a,k ' 'l a ,k ~ X a,k ^ 4 «a,fc • \ u L I 

The system of coupled differential equations can then be 
written in the form (cf. Eq. (A2) of |l6[) 



xf (*)=W*(t)X^(i) 



(/3)/ 



(D2) 



where 



X 



(0) 



£n h—X„ k v 0,k ■■■ V n Iaa . x .kyo,k ■■■2/n max .fc / 

(D3) 



The matrix W fc (t) is given by Eq. (A4) of g6| but 
here indices start at zero. The number of produced 
gravitons can be calculated direct ly from the solutions 
to this system using Eqs. (3.28) and ( 3.32| ). Note 
that for a given truncation parameter n max the above 
system of size 4(n max + I) x 4(n max + I) has to be 
solved n max + 1 - times, each time with different initial 
conditions ( |3.3S| ). 

The main difficulty in the numerical simulations is that 
most of the entries of the matrix Wfc(i) [Eq. (A4) of 
]f6| ] are n ot kn own analytically. This is due to the fact 
that Eq. ( 2.40 ) which determines the time-dependent 
KK masses mi(t) does not have an (exact) analytical 
solution. Only the 00-component of the coupling matrix 
M a p is known analytically. We therefore have to deter- 
mine the ti me-de pendent KK-spectrum {rrii(t)}"^ x by 
solving Eq. (|2.40 ) numerically. In addition, also the part 



M-j [Eq. (BE)] has to be calculated numerically since 



the integral over the particular combination of Bessel 
functions can not be found analytically. 

We numerically evaluate the KK-spectrum and the 
integral for discrete time-values ti and use s pline 



routines to assemble Wfc(i). The system (D2) can 
then be solved using standard routines. We chose the 
distribution of the ti's in a non-uniform way. A more 
dense mesh close to the bounce and a less dense mesh at 
early and late times. The independence of the numerical 
results on the distribution of the t$ 's is checked. In order 
to implement the bounce as realistic as possible, we do 
not spline the KK-spectrum very close to the bounce 
but re-calculate it numerically at every time t needed in 
the differential equation solver. This minimizes possible 
artificial effects caused by using a spline in the direct 
vicinity of the bounce. The same was done for MfJ but 
we found that splining My when propagating through 
the bounce does not affect the numerical results. 
Routines provided by the GNU Scientific Library (GSL) 
p9[ have been employed. Different routines for root 
finding and integration have been compared. The code 
has been parallelized (MPI) in order to deal with the 
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FIG. 33: Comparison of the final KK-graviton spectrum 
J\fnk • with the expression d n ,fc(iout) describing to wha t accu- 
racy the diagonal part of the Bogoliubov relation (D4) is sat- 
isfied. Left panel: y s = 3, k — 0.1, Vb = 0.03 and n max = 100 
[cf. Fig. ||. Right panel: y s = 3, k = 30, v b = 0.1 and 
n max = 100 [cf. Fig. M. 



intensive numerical computations. 

The accuracy of the numerical simulations can be 
assessed by checking the validity of the Bogoliubov 
relations 

E [^a,*(*)^,fc(*) " B2a,*(*)*W(*)] = ««7 (° 4 ) 

E l/WW^kW ~ BJa,*(*)-^/S7,*(*)] = 0- (D5) 

/9 



In the following we demonstrate the accuracy of the nu- 
merical simulations by considering the diagonal part of 
(D4). The deviation of the quantity 



da,k(t) = 1 - E [l-WWI 2 - I^WWH (D6) 



from zero gives a measure for the accuracy of the 
numerical result. We consider this quantity at final 
times tout and compare it with the corresponding final 
particle spectrum. In Fig. [33| we compare the final KK- 
graviton spectrum with the expression d n ,fc(iout) 
for two different cases. This shows that the accuracy 
of the numerical simulations is very good. Even if the 
expectation value for the particle number is only of order 



10~ 7 to 10~ 6 , the deviation of d n> k(t ut) from zero is at 
least one order of magnitude smaller. This demonstrates 
the reliability of our numerical simulations and that we 
can trust the numerical results presented in this work. 

APPENDIX E: DYNAMICAL CASIMIR EFFECT 
FOR A UNIFORM MOTION 



We consider a real massless scalar field on a time- 
dependent interval [0, y(t)]. The time evolution of its 
mode functions are described by a system of differential 
equations like ( 2.49 ) where the specific form of M a p de- 
pends on the particular boundary condition the field is 
subject to. In Jl5|, [TtJ a method has been introduced to 
study particle creation due to the motion of the boundary 
y(t) (i.e. the dynamical Casimir effect) fully numerically. 
We refer the reader to these publications for further de- 
tails. 

If the boundary undergoes a uniform motion y(t) = 1+vt 
(in units of some reference length) it was shown in j57], pgfl 
that the total number of created scalar particles diverges, 
caused by the discontinuities in the velocity at the begin- 
ning and the end of the motion. In particular, for Dirich- 
let boundary conditions (no zero mode), it was found in 
Hi that (0,in|JV° ut |0,in) oc v 2 /n if n > 6 and v < 1. 
Thereby in- and out- vacuum states are defined like in 
the present work and the frequency of a mode function 
is given by ui n — irn , n = 1,2, ... . In Figure [m] we 
show spectra of created scalar particles obtained numer- 
ically with the method of JItJ for this particular case. 
One observes that, as for our bouncing motion, the con- 
vergence is very slow since the discontinuities in the ve- 
locity lead to the excitation of arbitrary high frequency 
modes. Nevertheless, it is evident from Fig. [54| that the 
numerically calculated spectra approach the analytical 
prediction. The linear motion discussed here and the 



brane-motion ( 2. IS ) are very similar with respect to the 
discontinuities in the velocity. In both cases, the total 
discontinuous change of the velocity is 2v and 2vb, re- 
spectively. The resulting divergence of the acceleration 
is responsible for the excitation and therefore creation of 
particles of all frequency modes. Consequently we ex- 
pect the same oc v 2 jui n behavior for the bouncing mo- 
tion (2.18). Indeed, comparing the convergence behavior 
of the final graviton spectrum for Vb = 0.01 shown in 
Fig. ^5| with the one of the scalar particle spectrum for 
v = 0.01 depicted in Fig. |34| shows that both are very 
similar. 
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